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CONVERSION  OF  KINETIC  ENERGY  OF  A  STATIONARY  PLASMA  STREAM 
INTO  ELECTRICITY  WITH  ISOTHERMAL  RETARDATION 
OF  THE  STREAM  IN  A  TRANSVERSE  MAGNETIC  FIELD 
V. M.  Saryche v 
(Novosibirsk) 

Recently  many  investigations  have  been  made  of  the  flow  of  plasma 
in  so-called  magnetohydrodynamic  generators.  In  Reference  [1]  was  con¬ 
sidered  the  steady  motion  of  a  one -dimensional  stream  of  an  incompres¬ 
sible  nonviscous  liquid  through  a  homogeneous  magnetic  field  in  a  chan 
nel  of  arbitrary  cross  section,  and  that  of  a  compressible  nonviscous 
liquid  in  a  channel  of  constant  cross  section  at  maximum  power  deliv¬ 
ered  to  the  external  circuit  per  unit  channel  length.  In  Reference  [2] 
there  was  considered  steady  motion  of  a  one -dimensional  stream  of  a 
compressible  nonviscous  liquid  in  a  channel  of  arbitrary  cross  section 
through  a  transverse  homogeneous  magnetic  field  with  the  velocity  of 
the  liquid  stream  and  the  intensity  of  the  electric  field  constant 
along  the  channel. 

Since  the  maximum  temperature  that  modern  materials  can  withstand 
and  the  minimum  temperature  at  which  plasma  can  be  produced  with  con¬ 
ductivity  sufficient  for  effective  interaction  with  a  magnetic  field 
are  approximately  equal,  it  is  of  interest  to  consider  the  flow  of  a 
plasma  in  a  magnetohydrodynamic  generator  with  a  plasma  temperature 
constant  along  the  channel. 

In  the  case  of  isothermal  retardation  of  the  plasma,  one  can  use 
for  generation  of  electricity  only  the  kinetic  energy  of  the  plasma,  . 
the  ratio  of  which  to  the  total  energy  of  the  plasma  stream  (without 
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the  ionization  and  dissociation  energy)  is  in  the  case  under  considera¬ 
tion 


(T-imvar 
(t  — i).w»/ar  +  * 


Here  M  is  the  Mach  number  for  isothermal  flow.  When  M  =  1  the  ra¬ 
tio  tj  =  (y  —  1)/(3y  —  1)  for  a  monatomic  gas  amounts  to  'approximately 
17#. 

*  *  * 

1.  System  of  equations.  We  consider  isothermal  motion  of  a  one¬ 
dimensional  stream  of  compressible  nonviscous  liquid  with  finite  con¬ 
ductivity  in  a  planar  channel  of  arbitrary  rectangular  cross  section 
through  a  transverse  magnetic  field.  The  channel  walls  parallel  to  the 
magnetic  field  are  electrodes  which  make  it  possible  for  a  potential 
difference  that  alternates  along  the  length  of  the  channel  to  be  formed 
on  them  ( sectionalized  electrodes,  separated  by  insulators).  The  elec¬ 
trodes  are  interconnected  through  an  external  load  in  such  a  way  that 
the  resistance  of  the  latter  is  arbitrarily  distributed  over  the  length 
of  the  channel,  while  the  current  in  the  plasma  flows  perpendicular 
to  the  stream  velocity  of  the  plasma  and  to  the  magnetic  field.  We  are 
considering  a  case  when  the  magnetic  field  induced  by  the  currents 
flowing  in  the  plasma  can  be  neglected  compared  with  the  external  mag¬ 
netic  field. 

Let  vf  =  (u,  0,  0),  7=  (0,  — j,  0),  if  =  (0,  0,  H).  The  motion  of 
the  plasma  can  be  described  in  this  case  by  the  following  system  of 
equations: 

pub  =  p0uQb0  (continuity  equation)  (1-1) 

p  =  pRTQ  (equation  of  state)  (1.2) 

+  (equation  of  motion  (1-3) 

*  projected  on  the 
channel  axis) 
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P»*  jjjf-  *■  — 7*  -f  ( energy  equation ) 


(1.4) 


— I  -)  (equation  of  the  elec-  (1.5) 

*  \  -hr* If*  ftlrftiHt^ 


*  '•  *'  trie  circuit) 

Here  r*  is  the  resistance  of  the  external  circuit  per  unit  area 
of  the  walls  (the  channel  electrodes).  We  shall  consider  two  depend¬ 


ences  of  the  plasma  conductivity  a  on  its  density  p 


(a)  o  =s  const,  (b)  e  — 


Vt 


(1.6) 


The  first  dependence  takes  place  when  the  plasma  is  made  of  a  gas 
with  high  ionization  potential,  with  small  additives  of  readily  ion¬ 
ized  matter,  and  the  plasma  temperature  and  composition  are  such  that 
the  main  gas  is  not  ionized  but  the  additives  are  almost  completely 
ionized  (3). 

The  second  dependence  takes  place  in  the  case  of  a  homogeneous 
weakly  ionized  plasma,  if  the  Interaction  between  the  electrons  and 
the  ions  can  be  neglected  compared  with  the  interaction  of  the  elec¬ 
trons  and  the  neutral  particles  (neSel  «  nnSen)* 

For  the  length  scale  we  use  the  width  of  the  channel  b^  at  the 
section  x  =s  0;  the  scales  for  the  pressure,  density,  and  velocity  of 
the  plasma,  density  of  the  current,  intensity  of  the  magnetic  field, 
conductivity  and  resistance  of  the  external  circuit  per  unit  area  of 
the  electrodes  we  use  the  values  of  the  following  quantities:  pQ,  pQ, 
uQ,  jQ,  HQ,  o0,  reQ  in  the  section  x  =  0.  We  shall  use  dimensionless 
variables  throughout  without  specially  designating  them. 

2.  Optimal  retardation  of  the  plasma  stream.  It  follows  from  (1.4) 
and  (1.5)  that  the  maximum  retardation  of  the  plasma  and  the  maximum 
energy  delivery  to  the  external  circuit  per  unit  length  of  the  channel 


occur  [1]  when 


re  -  b/a 


(2.1) 
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Unlike  ordinary  generators,  when  the  energy  released  in  the  in¬ 
ternal  resistances  is  parasitic,  in  a  magnetohydrodynamic  generator 
this  energy  remains  in  the  plasma  and  can  be  delivered  to  an  external 
circuit.  (In  the  case  of  isothermal  retardation,  the  energy  delivered 
to  the  external  circuit  depends  only  on  the  values  of  the  initial  and 
final  velocities  and  does  not  depend  on  the  distribution  and  value  of 
re.)  Therefore  the  mode  (2.1)  is  very  favorable  for  a  magnetohydrody¬ 
namic  generator. 

(A)  We  consider  a  case  when  H  =  1  and  a  =  1.  In  this  case  the  sys 
tem  of  equations  (l.l)-(1.6)  is  closed.  Solving  it,  we  obtain 


b  —  rm 

db 

dx 


M 

'=  AT.* 

M, 

=  -g-  exp  - 

<?,(*/»  +  !) 

2  M* 


P=  P 
M*  —  M* 


“  exp  [ 


M*  —  A/„» 


w. 

--qt)cx  p 


]  (*-t) 

a n-Mf 


<7o 

wv 


dM 


(2.2) 


p* 


Here  M  is  the  Mach  number  for  the  isothermal  flow  and  Aq  and  Qq 
are  arbitrary  dimensionless  parameters. 

In  the  case  when  a  =  1/ the  solution  has  the  form 

b  = 


p  =  p  =  exp 


M* —  A#0* 


M,  M*-M* 


M 


exp 


M 

I  “  Si  exP 


Mt  ™o- 

r*=  AT  exP  — 


2 

MJ  —  M* 


(2.3) 


x  exp  [4-  (M*  -  A/.*)]  dAf 

5- = «P  [4*  w  - 


(B)  Let  us  consider  the  problem  of  the  maximum  isothermal  re  tarda 
tion  of  the  plasma  stream  with  the  potential  difference  on  the  channel 
wall-electrodes  remaining  constant  along  the  channel 


Jre  =  1 

For  o=l  the  solution  has  tie  form 


(2.4) 


M  If* —  AfJ  /  1 

/=TSrexp — 2 —  (*-t cxp  — - 


■  M* 


dM 


(2.5) 
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p  =  p  mi  H  =  ^  =  e  xp  -------  ,  p  =  pt=//-*e*p-^l  as  At-** 

The  angle  of  channel  expansion  will  in  this  case  be  determined  by 
the  formula 


M  _  Q,  M*  +  1 

"S'  =  1 - w~ 


(*£",-*  as  *-«) 


For  a  =  l/<y/p”  the  solution  has  the  form 


^  =  p  =  //  =c  =  rxp 


/  =  jmT”P - ? -  (*-t) 

Af. _ __  2Af,  f* _ 


(2.6) 


(2.7) 


,  .X  = 


IS")  exP  T" 


«W 


The  angle  of  channel  expansion  is  determined  in  this  case  by  the 
formula 

m  q,  w*+i _  «!-**  /*  _ _  „\ 

TF  =  ”2 - Hr— exP  l  \~dT~~  as  M~0)  {2.8) 


3.  Retardation  of  the  plasma  In  a  channel  with  constant  expansion 
angle. 

(A)  Let  H  =  1  and  0  =  1.  We  then  get  from  (1.5) 

r*“-sfc-(TO~l)(1+x*)  (*-2t«-r)  (3.1) 

Here  8  Is  the  angle  of  channel  expansion.  For  J  we  obtain  In  this 

case 

1 "  i\U.  jpcrt  (- 1  +  [l  +  4%-fc  (3.2) 

When  M  -*•  1  we  have  J  -*  <*>;  it  is  therefore  Impossible  to  makse  a 
continuous  transition  through  the  point  M  .=  1. 

When  M  <  [1  —  Qq(1  +  xx)2 311(1  when  Aq(M  —  1)  <  0  the 
problem  considered  has  no  solution.  When  x  =  0  we  have 

<3-3> 

It  follows  therefore  that  the  problem  has  no  solution  when 

x/l  0*1X1  A,|<max  {1.  x/|0*|> 

Solving  the  Initial  system  of  equations  with  respect  to  M,  we  ob- 


-5 


tain 


i  dM 
ST  dx 


Q.<\+Xr)U* 
ZA-Ui. 


+  Xr)  ft  . 

it(W*'=TF\  ^ 


(J-  Xx)  (.»/•-«) 

+i«+*-**'san 

It  follows  therefore  that  dM/dx  <  0  only  If 

^  Vx+ns 
M<  iTv 

In  the  case  when  a  =  1 /yfp  we  have  for  J  and  M 

/  -staf  *&£  {-"  [>+<«-&*.•(£■)'  ^r]"} 


(3.4) 


1  4M 


_ X  0.  /  W  \V.  (I  +  X«)%  li  J. 

If  ds  ~  («•— 1)<1  +  X*)  ■*  2/4,  \1?W  (W»  — 1)*  \ 


J3-5) 


(B)  Let  us  consider  the  case  when  Jr  =  1  and  a  =  1,  Solving  the 

w 

system  of  equations  (l.l)-(1.6)  and  (2.4)  we  obtain 
1  +  x*  =  exp  5  A/*  {  - 1  +[  1  +  ^%rr*: 


P  =  P 


w. 


„  _  r,  _ 1 _ i _ 1  (i  -  a.)  u, 

’  "“L  M'-\  M  {i+x*)4M  ldz\  <l+x*)* 


W  (1+X*) 


/LL  _1  +  <v«\ 

’~r,  -/«,[*»,  t  +  X*  J  U.-'+’V' 


(3.6) 


The  problem  has  a  solution  when  0  <  A~  <  1. 

Isothermal  retardation  of  the  plasma  in  channels  of  constant 
cross  section  and  in  narrowing  channels  (x  <  0)  is  impossible. 

4.  Retardation  of  a  plasma  in  a  homogeneous  magnetic  field  with 
the  potential  difference  on  the  wall-electrodes  of  the  channels  con¬ 
stant  along  the  length  of  the  channel.  For  a  =  1  the  solution  of  this 
problem  has  the  form 


rid' 

1  /.,.p 

0- 

J# 

U-(t-A,)pi  | 

[*- 

i  + 


N, 

*.  c 


exp  -Llf5£  *,/ 


(4.1) 


A  o-^)1  ftwpt  V.  Af*)|  (1  +  /4,  (A/1  - 1) ««p  IV, (Uf  -  M*)D 

**  “  (I  +  /4,«xp  (*/,  (M»-  **)))» 
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The  problem  has  a  solution  when  0  <  AQ  <  1.  If  a  =  l/y^p",  then 
the  solution  has  the  form 

' “ [* —  '»•*•*►  (— -)]  nr*-  *“  Sr 


>~T' 


EVTS-"-"-^'1  (-4- *+■?*•) 

M9  KM  It  •  -  u  .  att 


j ,  *>  r _ 


exp 


(4.2) 

0 


Qb(l  —  tOv*  J  r  2 
and  exists  also  only  In  the  range  0  <  Aft  <  1. 
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GENERATION  OF  ELECTRIC  CURRENT  BY  A  PLASMA  FLOWING  IN  A  TUBE 

E. I.  Andriankin 
(Moscow) 

The  simplest  generator  scheme  in  the  form  of  a  tube  of  rectangular 
cross  section  was  first  investigated  theoretically  in  [1],  where  the 
efficiency  of  such  a  generator  was  estimated  in  the  case  of  flow  with 
constant  conductivity. 

It  is  known  that  the  conductivity  of  an  ionized  gas  depends  ex¬ 
ponentially  on  the  temperature,  which  varies  noticeably  along  the  tube. 
Therefore  for  a  quantitative  estimate  of  the  efficiency  it  is  neces¬ 
sary  to  take  into  account  the  temperature  dependence  of  the  conductivity. 

In  the  present  investigation  we  estimated  the  electric  efficiency 
of  a  generator  with  conductivity  corresponding  to  a  mixture  of  nitrogen 
and  1 %  potassium  vapor.  In  the  calculation  we  used  the  interpolation 
formulas  obtained  by  G. S.  Aravin  and  V. P.  Shevelev  for  the  dependence 
of  the  conductivity  on  the  pressure  and  the  temperature  [2],  The  re¬ 
sults  of  the  calculations  are  in  qualitative  agreement  with  the  data 
of  [3 ],  where  the  flow  of  cesium  vapor  was  investigated  in  a  magneto¬ 
hydrodynamic  current  generator.  The  presence  of  friction  against  the 
tube  walls  influences  in  some  cases  noticeably  the  flow  of  the  plasma 
and  will  also  be  accounted  for  in  the  analysis  of  the  efficiency. 

In  solving  the  problem  we  disregarded  several  factors:  the  aniso¬ 
tropy  of  the  conductivity  [7,  8],  the  terminal  losses  at  the  closed 
points,  the  heat  dissipation  of  the  walls,  and  finally,  the  losses  to 
electric  resistance  in  the  near-electrode  boundary  layer.  Therefore 
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the  true  values  of  the  efficiency  are  lower  for  the  examples  considered 
In  the  present  paper.  We  note,  however,  that  to  find  the  most  suitable 
generator  mode  there  are  still  several  possibilities  left:  variation 
of  the  external  magnetic  field  and  variation  of  the  cross  section  of 
the  tube  along  the  length  of  the  channel  [9],  use  of  sectlonallzed 
electrodes  to  pick  off  the  Hall  current,  etc. 

#  #  * 


In  the  case  of  isotropic  conductivity  [4],  the  flow  of  a  plasma 
in  a  current  generator  (Fig.  1)  whose  horizontal  walls  are  maintained 
at  a  constant  voltage,  and  the  magnetic  permeability  and  the  dielectric 
constant  of  the  plasma  are  close  to  unity,  is  described  by  the  system 
of  equations 


■5r  +  TT“0'  p* ■£-  +  p*’-£-+-£-  +  “¥L  +  piS 
P“  +  P*  %  +  "Sr—  -p~  =  °t  »  -  »  /*'■<*>- 

it  [  p“  [^i~  +  w)  -  -w *•*•] + 

+  [  Pp  +  w)  +  TT £*<f«]  ■ 


Fig.  1 


Here  u  and  v  are  the  velocity  com¬ 
ponents,  2  the  pressure,  p  the  density, 
r(x)  the  external  resistance,  and  dF^/dx 
denotes  the  friction  force,  which- depends 
on  the  value  of  the  magnetic  field,  and 
which  Is  averaged  over  the  cross  section. 


The  boundary  conditions  for  Eqs.  (1)  will  be  the  equations  of  con¬ 
tinuity  on  the  planes  ±Vq/2  for  the  components  of  the  magnetic  field, 
ahd  the  vanishing  of  the  normal  component  of  the  stream  velocity.  The 


initial  parameters  of  the  stream  (uq,  Pq,  p0)  at  the  section  x  *  0  can 
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be  specified  arbitrarily  only  in  such  a  range  of  variation  for  which 
the  flow  out  of  the  tube  remains  in  the  sonic  mode.  In  the  case  of  sub 
sonic  flow,  the  parameters  at  the  entrance  to  the  tube  are  connected 
with  its  length  1  and  with  the  external  pressure  p#.  This  connection 
is  obtained  by  solving  the  boundary  problem.  To  estimate  the  various 
terms  in  Eqs.  (l)  we  seek  the  solution  in  series  form 

u  =  i»,  (*)  +  !»,(*,  y)  + ... ,  p  =  p,  (x)  +  pt  (*.  y)  +  •  •  • . 

p  =  r,(x,  >/)  +  ...,  P  =  Pi(*)  +  p»(x,  y)  +  ...,  (2) 

H y  ~  //*  +  lift  (x,  y) *f  •  •  •  i  Hx  —  f?tx  (*i  •• 

After  solving  the  equations  for  u^,  p^,  and  p ^  we  determine  the 
current 

From  the  current  distribution  we  can  calculate  the  vortical  mag¬ 
netic  field 

H  =  -rS^rfw*  R  =  (*-*')  i«  +  Or  - </’) •,  +  (*- z') h. =  dx’dy'dz’  (3) 

Volume  integration  shows  that  Hgx  has  a  maximum  in  the  case  of 
flow  between  parallel  planes  (Zq  -*•  <*>) 

(4) 

Taking  (4)  into  account  and  equating  the  inertial  terms  in  the 
equations  of  motion  to  the  Lorentz  forces,  we  note  that 

Pi  ~  Pi«i*fl* 

It  then  follows  from  the  continuity  equation  that  Vg  <  u^q/1.  If 
Yq/1  «  1,  then  the  velocity  component  Vg  can  be  neglected  in  (1)  even 
if  account  is  taken  of  the  influence  of  the  induced  field  on  the  main 
stream.  Practical  interest  attaches  to  the  case  when  the  magnetic  Rey¬ 
nolds  number  R#  is  small.  We  \ntroduce  the  dimensionless  quantities 
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.  _  ///.*=. 

* 


X  = 


'  -  *•/ 


Ti%  # 


•  9* 


(5) 


•»«-=■ 


Omitting  terms  proportional  to  R#,  we  write  Eqs.  (1)  in  dimension 
less  form 

+  +a,,l'~ x,  +  0  £  =  0'  ®*  =  ®» (6) 

■sr[x+<T^i?]+,w'|K““,  =  0’  Pr=QW‘  QV=i 

The  parameter  tj  represents  the  ratio  of  the  ponderomotive  force 

A 

to  the  inertial  force  Pquq  /*J  the  parameter  k  characterizes  the  ra¬ 
tio  of  the  voltage  drop  on  the  external  circuit  to  the  electromotive 
force  in  the  considered  section  of  the  channel,  and  ranges  from  zero 
to  unity,  when  r  changes  from  zero  to  infinity. 

The  boundary  conditions  for  (6)  are  written  in  the  form 

P  =  1,  V  =  1  for  tj  =t  0 
P  .  P,  f or  tj  =  6,  ji#  <  1  (7) 

P  =  for  tj  =  6,  ji#  »  1 

(the  asterisk  pertains  to  the  quantities  at  the  output  from  the  tube). 
Analyzing  Eqs.  (6),  we  can  establish  that  when  ji  <  1 

TT>0’  £>°.  £<0.  £<0.  £<0 

In  the  case  when  ji  >  1  the  signs  of  the  inequalities  are  reversed 
Let  us  estimate  now  the  maximum  possible  generator  efficiency  as¬ 
suming  that  there  is  no  viscosity  and  the  magnetic  Reynolds  number  is 
small.  In  this  case  Eqs.  (6)  admit  of  the  Important  Integral 

•  c.^T+»..v-.(^+x)  (8) 

Relations  (8)  are  sufficient  for  an  estimate  of  the  efficiency, 
which  we  take  to  mean  the  ratio  of  the  difference  in  the  energy  fluxes 
at  the  input  and  output  of  the  tube  to  the  energy  flux  at  the  input. 
The  expression  for  the  efficiency  can  be  written  in  the  form 
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Therefore  the  efficiency  tends  to  its  maximum  value  when  is 
minimal,  i.e.,  when  k  -*  1.  The  dependence  ( |jl ,  l)  for  y  —  1.4  is 

shown  in  Fig.  2  (curve  a-b). 

p 

If  the  external  pressure  P#  turns  out  to  be  larger  than  p.  V^, 
then  the  outflow  mode  will  be  subsonic.  For  each  0  <  P#  <  1  there  ex¬ 
ists  a  p.Q  <  1  such  that  when  p.  <  pQ  the  outflow  mode  is  subsonic,  and 
when  p.  >  p.Q  the  outflow  mode  remains  sonic  even  when  /c  =  1.  From  the 

p 

condition  p.Q  =  P*  we  determine 

a y  «* —  T*  *  t*  (t*  ~  i)  ,  . 

K.’=  -  ~ L-  «  - f. (T* - 1)  +  1  ( 12) 

The  velocity  at  the  output  of  the  tube  in  the  case  when  p.  <  p.Q  is 
given  by  the  formula 


Since 


(13) 


we  have 


as  »**i,  i»-»o 


(14) 


_  rr-Dd-y.) 

▼max  7 

Let  us  prove  that  the  lines  9(1,  p.,  P*)  lie  when  p.  <  p.Q  below  the 
curve  (a-b)  on  Fig.  2.  Assume  that  there  exists  a  region  where  <p  >  9^. 
By  virtue  of  the  conditions 

f  (».  m*.  p.)  -=  f»(i.  M.  f  (i, »,  P ,)  <  *  (l,  0) 

there  exists  an  intersection  point  0  <  p.g  <  p.Q,  where 
f  (*.  I»t.  P»)  =  f»0.  M»)  *=  f  (1,  M«.  **).  Pt  “  |k) 


Inasmuch  as  dPg/3p.g  >  0  we  have  Pg  <  P#.  On  the  other  hand,  the 
derivative  69/dP  <  0  is  a  monotonic  function  and  consequently  the  effi¬ 
ciencies  cannot  be  equal  at  two  different  points  with  respect  to  P  for 
one  and  the  same  p.  =  p.g,  which  proves  our  statement.  The  curves  9(1,  p., 
P#)  for  the  cases  P#  =  0.1,  0.3,  and  0.5  are  shown  in  Fig.  2. 

Let  us  consider  now  the  conditions  for  the  maximum  efficiency  in 
the  case  of  supersonic  velocity  at  the  inlet  to  the  tube,  when  two  out¬ 
flow  modes  are  possible:  sonic  and  supersonic.  An  analysis  of  Formula 
(il)  Bhows  that  there  exists  a  number 


_ -«l+  V ■*?  /«i-(T-D(T|»*  +  1)  \  /1C\ 

^  * -  L-<t+i>(2-t>W  K0) 

for  which  >  *  if  k  <  k+  and  =  <c+  when  k  =  k+.  From  physical  con¬ 
siderations  V  >  k.  Therefore  for  k  >  »c+  only  supersonic  outflow  is  pos¬ 
sible  (the  minimum  velocity  at  the  output  is  k). 

For  the  sonic  outflow  mode  dV^/dic  <  0,  therefore  the  maximum  value 
of  9^  la  attained  when  =  ic+.  The  following  inequality  holds  true 


/**\  _  x  |  P  (V)  -  n^i 

,  ~  tTi  T  * 77 (t  <  0  <*•>*>  ' 

Therefore  to  determine  the  maximum  efficiency  in  the  case  of 
supersonic  outflow  it  is  necessary  to  put  in  Formulas  (8)  and  (9)  V 
and  obtain  the  value  of  jcq  which  gives  the  extremum  of  9 
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■1  ±  TX* 


(T+l)|»* 

Comparing  (15)  with  (16)  we  obtain 


,  »i-%- (t - ») [ 2(r,-ir(iV:f/;h -*)?>? ~ *]  (l6) 


x.-x+=^|m«(T-1)(2-T*H  (2-r)(2|**-l)|  .  M>0>  (17) 

It  follows  from  (17)  that  kq  >  *c+  when  7  =  1.4,  and  therefore  the 
value  of  which  leads  to  a  maximum  efficiency  equal  to  <p(«c0)  lies  in 
the  supersonic  region.  Figure  2  shows  the  curve  OC  calculated  from 
Formula  (16)  for  the  case  y  =  1.4  (the  scale  is  indicated  on  top  of 
the  plot). 

Let  us  estimate  the  influence  of  the  friction.  Following  [5],  we 
assume  that  the  friction  force  averaged  over  the  tube  cross  section  is 
proportional  to  the  velocity  head,  i.e., 

dF_  _  .  «■ 

dx  "  *  2» 

The  coefficient  £,  which  depends  on  the  magnetic  field  [6],  on 
the  turbulence  of  the  stream,  and  on  the  heat  conduction  on  the  wall, 
will  be  assumed  to  be  a  constant  quantity  obtained  by  experiment. 

Eliminating  the  variable  q  from  (6),  we  arrive  at  an  equation  re¬ 
lating  the  pressure  with  the  velocity 

dP  ...  («»V  —  T3,x)J*  +  (t  —  — 

~dV --V-  e,l'a  +  |'8|K(l_2T)  +  ^8|(T_|) 


djL 

dV 


TO!  (V  —  x)  +  al 


wy  (18) 


If  o1  is  constant,  then  Eqs.  (18)  can  be  integrated  in  quadra¬ 
tures.  In  this  case,  however,  it  is  possible  to  obtain  a  simple  expres¬ 
sion  for  the  efficiency  only  when  a  is  small.  Urns,  for  example,  for 
p.  -*■  0,  k.  =  1,  using  the  perturbation  method,  we  obtain 

v(0,l)  =  IriL[  1  —  /*,  +  a  In  _^  + !^!<1  _*>,)«]  (19) 

It  follows  from  (19)  that  the  friction  must  be  taken  into  account 
even  for  small  a,  since  its  influence  on  the  efficiency  is  proportional 
to  a  In  a.  It  is  possible  to  estimate  the  Influence  of  the  viscosity 
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Pig.  3 


more  rigorously  using  as  an  example  the  flow  of 
an  Incompressible  plasma  between  parallel  planes 
(Zq  -*■  <*>)  in  an  external  magnetic  field  [4]  (the 
Hartmann  problem).  The  expression  for  the  effi¬ 
ciency  in  the  maximum-power  mode,  which  in 
this  case  must  be  taken  to  mean  the  ratio  of 
the  work  done  in  the  external  circuit 

I  Mft 

J  \  r  {x)j&ydz 
to  the  pushing  work 

a  t,  x 

4* 


has  the  form 


_  cch«  —  she 
2(ech«+*h®)  * 

0  =  0, 


(20) 


Here  X  is  the  viscosity  coefficient. 

In  order  to  calculate  the  maximum  efficiency  for  a  given  generator 
length  (6  =  const)  it  is  necessary  to  Integrate  Eqs.  (18)  numerically. 
It  was  assumed  in  the  calculations  that  the  plasma  conductivity  corres¬ 
ponds  to  thermal  ionization  of  a  mixture  of  nitrogen  with  l£  potassium 
vapor  and  is  represented,  in  accordance  with  the  calculations  of  [2], 
in  the  form 


°"=0*  Y~^ cxP~ri='^' exP [*.(  (21) 

“  «•  *'*»  *  *  y- ,  Pi  —  1  atm 

«  =  2.6-lO*cK,  o,  =  !).05*1(F  l/ohm—  cm. 


For  subsonic  entrance  velocities,  the  parameter  b  was  assumed 
equal  to  8.5,  corresponding  to  an  initial  temperature  in  the  plasma  of 
approximately  3000°K. 
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Figure  3  shows  the  dependence  of  the  efficiency  and  of  the  dimen¬ 
sionless  tube  length  6,  for  which  the  sonic  outflow  mode  is  attained, 
on  the  parameter  c  at  o  =  0  and  for  p  varying  from  0.1  to  0.6.  As  k  -*  0 
the  curves  6  =  6(k,  p)  tend  to  a  finite  limit,  and  as  k  -*  1  they  in¬ 
crease  without  limit,  since  the  equations  (18)  have  a  singularity  in 
the  case  a  =  0,  k  =  1. 


Fig.  4 


Fig.  5 


Fig.  7 


Figure  4  shows  an  analogous  dependence  for  the  case  a  —  0.1.  When 
a/0  Eqs.  (18)  do  not  have  any  singularities  and  the  6(k,  p)  curves 
tend  as  /c  -*■  1  to  a  finite  limit.  This  is  quite  obvious,  for  owing  to 
friction  (even  in  the  case  of  flow  without  a  magnetic  field)  the  sonic 
mode  will  be  attained  in  a  tube  of  finite  length  [5]« 

Figure  5  shows  the  dependence  of  the  efficiency  on  the  parameter 
jc  for  the  following  cases,  respectively:  (a)  a  =  0,  p  =  0.3;  (b)  a  = 

=  0.1,  p  =  0.3*  and  (c)  a  =  0.1,  p  =  0.15;  the  values  of  the  parameter 
6  are  indicated  on  the  curve.  From  a  comparison  of  Figs.  4  and  5  it 
follows  that  the  efficiency  for  constant  p  increases  with  increasing  6, 
but  always  remains  smaller  than  in  Fig.  2;  the  maximum  value  of  the 
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efficiency  for  constant  p  and  for  constant  6  is  attained  for  suffi¬ 
ciently  large  6  in  the  case  of  a  sonic  mode  at  the  tube  output  (curves 
OA  on  Pig.  5).  For  a  tube  with  specified  value  of  the  parameter  6  = 

=  const,  there  exists  a  fully  defined  Mach  number  at  the  input  to  the 
tube,  at  which  the  efficiency  is  a  maximum.  This  statement  is  illus¬ 
trated  in  Fig.  6,  which  shows  the  dependence  of  the  maximum  efficiency 
on  the  number  p  for  the  case  6  =  10  and  a  —  6.1. 

Figure  7  shows  the  distribution  of  the  plasma  parameters  along 
the  length  of  the  tube  x^  =  x/l  for  the  case  p  =  0.2,  as  0.1,  k  «  1. 
It  is  seen  from  the  plot  that  the  conductivity  changes  by  one  order  of 
magnitude  along  the  tube.  In  order  to  compare  the  efficiency  of  the 
generator  with  subsonic  and  supersonic  velocity  at  the  input  to  the 
tube.  Figs.  8  and  9  show  the  dependence  of  qp  and  on  k  for  the  cases 
1.5  <.  a  =  0  and  also  the  dependence  of  the  efficiency  on  k  for 

p^  =  2  and  different  values  of  6g.  It  was  assumed  in  the  calculation 
that  the  slowing -down  temperature  at  the  inlet  of  the  tube  is  approxi¬ 
mately  3000°K,  and  therefore  the  initial  values  of  OgQ  turn  out  to  de¬ 
pend  on  the  Mach  number 


**  "  (^T1  W)  . 


(22) 


The  dashed  curves  on  Fig.  8  correspond  to  values  k  -  *+.  To  the 
left  of  k  -  k+  the  ourvtB  tend  to  a  finite  limit.  Between  the 

dashed  ourve  and  the  curve  a-b,  corresponding  to  Eq.  (16),  we  have 
&2  "•  However,  if  friction  Is  taken  Into  account,  then  will  al¬ 

ways  be  finite.  In  the  case  of  a  supersonic  flow  out  of  the  tube,  a 
complicated  system  of  shocks  arises.  We  note  that  the  pressure  which 
would  occur  In  the  case  of  slowing  down  of  the  escaping  stream  In  a 
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direct  shock  wave  should  be  larger  than  the  external  pressure.  If  this 
is  not  satisfied,  then  a  shock  wave  is  produced  in  the  middle  of  the 
tube  in  such  a  place  that  the  internal  and  external  pressures  become 
equal  at  the  output  of  the  tube. 
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Fig.  8  Pig.  9 

it  follows  from  Fig.  8  that  for  large  jx  and  for  6  -*■  «.  the  effi¬ 
ciency  of  a  supersonic  generator  becomes  higher  than  in  the  case  when 
ix  <  1.  On  the  other  hand,  if  we  compare  generators  of  equal  length  for 
the  same  slowing-down  parameters,  then  we  have  from  (22)  <5g  «  6; 
therefore  in  the  cases  of  practical  interest  the  subsonic  modes  are 


more  advantageous  than  the  supersonic  ones.  Calculations  have  con¬ 
firmed  the  need  for  taking  into  account  the  dependence  of  the  conduc¬ 
tivity  on  the  temperature  and  on  the  pressure.  Thus,  for  the  case  6  = 
_  io,  n  =  0.3,  a  =  0.1  the  maximum  of  the  electric  efficiency  is  «10#, 
and  in  the  case  of  flow  with  constant  conductivity  without  friction 
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THERMAL  IONIZATION  AND  ELECTRIC  CONDUCTIVITY 
OP  SOME  COMBUSTION  MIXTURES  AND  PRODUCTS 

G.S.  Aravin,  V. P.  Shevelev 
(Moscow) 

On  the  basis  of  the  thermal  nature  of  Ionization  (the  Saha  equa¬ 
tion)  ,  a  calculation  Is  made  of  the  electron  concentration  and  elec¬ 
tric  conductivity  of  various  inert  gases  (argon,  nitrogen,  helium) 
with  additives  of  the  most  easily  ionized  substances  (cesium,  potas¬ 
sium,  sodium)  in  the  temperature  range  1500-4000°K,  and  also  of  the 
combustion  products  of  some  gaseous  highly  calorific  fuels  containing 
the  same  additives  at  atmospheric  pressure.  The  dependence  of  the  elec¬ 
tric  conductivity  of  the  mixture  on  the  temperature,  pressure,  type  of 
gas,  and  amount  of  additive  is  analyzed.  The  calculation  is  made  with¬ 
out  account  of  the  influence  of  the  electric  and  magnetic  fields  on 
the  concentration  of  the  electrons  and  on  the  electric  conductivity. 

Hie  calculation  data  must  be  regarded  as  approximate. 

*  #  * 

1.  Thermal  ionization  of  a  mixture  of  gases.  At  sufficiently  high 
temperatures,  at  which  dissociation  and  ionization  of  the  mixture  com¬ 
ponents  takes  place,  an  exact  determination  of  the  total  composition 
of  neutral  and  charged  mixture  components  necessitates  simultaneous 
solution  of  the  system  of  chemical  and  ionization  equilibrium  equa¬ 
tions.  However,  as  shown  by  preliminary  calculations,  at  sufficiently 
low  degree  of  ionization  it  is  permissible  to  neglect  the  latter  and 
to  solve  first  the  chemical  equilibrium  system  and  then,  vising  this  as 
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the  final  chemical  composition  of  the  neutral  particles,  the  equillb- 
rium  ionization  and  ion-production  processes  are  considered  independ¬ 
ently.  We  shall  be  interested  henceforth  in  the  region  of  temperatures 
and  pressures  where  the  mixture  components  are  essentially  singly  ion¬ 
ized,  and  there  is  practically  no  secondary  ionization. 

A  real  mixture  consists  in  the  general  case  of  several  components 
capable,  on  the  one  hand,  of  being  ionized  and,  on  the  other  hand,  of 
forming  negative  ions;  therefore  in  the  derivation  of  the  general  re¬ 
lations  for  the  thermodynamic-equilibrium  electron  concentration  we 
shall  take  into  consideration,  unlike  the  preceding  investigations 
[1-6],  loth  ionization  and  production  of  negative  ions. 

Assume  that  under  constant-pressure  conditions  the  equilibrium 
ionization  and  i  cnbination  of  each  component  A  of  the  type  1,  com¬ 
prising  the  mixture,  occurs  in  accordance  with  the  relation  A^  A^+  + 
4  e,  while  the  process  of  formation  and  decomposition  of  the  negative 
ions  B"  of  type  j  obeys  the  relation  Bj~  Bj  4-  e.  Then,  expressing 
the  equilibrium  constants  in  terms  of  the  partial  pressures  of  the 
corresponding  reagents,  we  obtain 


A'«  = 


Ri 


PjP- 

pf 


(l.l) 


Prom  the  matter  and  charge  conservation  laws  and  from  the  equality 
of  the  total  pressure  of  the  ionized  mixture  to  the  sum  of  the  partial 
pressures  of  the  components  we  get 

*•= 5- (*+**)•  pf-Zu>,+pn  (i-2) 

2p*+  “P.+Sp*".  p  =  2(p.+ pO  +  2 (Pi  +  pf) + p.  ( 3) 

Here  p®  and  Pj  are  the  initial  (prior  to  the  ionization  and  forma¬ 
tion  of  the  negative  ions)  partial  pressures  of  the  ionized  and  elec¬ 
trically  negative  components  of  the  mixture ,  mQ  and  m  are  the  initial 

and  final  number  of  moles  of  mixture,  and  p  is  the  partial  pressure 

© 
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TABLE  1 

Values  of  the  Ionization  Potentials  and  the 

Ratio  of  the  Statistical  Weights  of  the  Qround 
State  2g+/g  for  Different  Components  (the  Litera¬ 
ture  Sources  are  Indicated  In  the  Square  Brackets) 


1 

1  j 

j  •’*  j  K 

Is* 

NO 

j  O. 

j  H/* 

OH 

II 

vi .»  3 

UJt 

3. OKI 
l»l 

1 

(19| 
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|I3| 

1 

n»i 
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|13| 
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13.505 

•l«l 

I* 

l**l 

1 

0 

CO. 

CO 

H, 

Ns 

EH 

■Ha 

-  '  >• 

3 

U*lt 

-o. 

_ nj _ 

13.700 

1151 

4 

I  Op«- 

ran- 

'■asm. 
|FVWP  ■ 

© 

5s 

15.427 

l*’l 

4 

|20J 

15.000 

(l«) 

4  . 

15.755 
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1)  Component;  2)  tentative;  3)  ev. 


of  the  electrons. 

Solution  of  this  system  of  equations  for  the  partial  pressure  of 
the  electrons  yields  the  following  expression 


-  /v  v  rfa  \ 


(1.*) 


Neglecting  here  the  value  of  pg  compared  with  £  (pe  «  p)  we  have 

U~5) 


We  note  that  this  assumption  Is  equivalent  to  considering  the  In¬ 
dicated  processes  at  constant  volume. 

Calculation  has  shown  that  In  the  case  of  high  Ionization  poten¬ 
tials  of  the  mixture  (for  example ,  in  combustion  products  without  ad¬ 
ditives  of  easily  ionized  substances.  Tables  1  and  2)  the  equilibrium 
constants  of  the  ionization  and  recombination  processes  are  many 
orders  of  magnitude  smaller  than  the  partial  pressure  of  the  electrons 

p  ,  and  the  latter  Is  in  turn  many  orders  of  magnitude  smaller  than 
© 
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the  equilibrium  constants  Kj  of  the  formation  and  decomposition  of  the 
negative  ions,  owing  to  the  small  values  of  the  electron-affinity  en¬ 
ergy  (0.87-1.735  ev).  Taking  this  into  account,  we  obtain  a  relatively 
simple  relation  for  the  determination  of  the  partial  pressure  of  the 
electrons  at  a  low  degree  of  ionization 

A-[2*«*V(l  +  2j£)]*  (1.6) 

which  in  the  case  when  no  negative  ions  are  produced  assumes  the  well- 
known  form  [4] 

(1*7) 

The  necessary  equilibrium  constants  contained  in  the  formulas  can 
be  calculated,  as  is  well  known,  from  the  Saha  equation,  which  for  a 
monatomic  gas  at  temperatures  at  which,  on  the  one  hand,  the  electron 
excitation  can  be  neglected  relative  to  the  ground  level*  and,  on  the 
other  hand,  the  fine  structure  of  the  ground  state  can  be  disregarded, 
has  the  form 

]g/r  =  _  + 

+  2.5  IgT-  6.48  +  lg^-\  (1*8^ 

ft 

Depending  on  the  process  considered,  K  denotes  here  the  equilib¬ 
rium  constant  of  the  ionization  and  recombination  processes,  K^,  or 
the  constant  for  the  formation  and  disintegration  of  negative  lonB,  Kj 
[atm],  I  is  the  ionization  potential**  or  the  electron-affinity  en¬ 
ergy  Ej  [ev],  while  g1  and  gg  are  the  statistical  weights  of  the  ground 
states  of  the  positive  ion  and  the  neutral  component,  or  the  electro¬ 
negative  component  and  the  negative  ion;  T  is  the  temperature  [°K]. 

As  applied  to  molecular  components,  Eq.  (1*8)  must  be  regarded  as 
approximate,  for  it  does  not  take  into  account  the  contributions  of 
the  rotational  and  vibrational  degrees  of  freedom  of  the  molecules. 
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TABLE  2 

Chemical  Composition  and  Temperature  T#  of  the  Combustion  Products  In  °K  of  Different 
Combustible  Mixtures  at  Pressure  p  =  1  atm  and  Initial  Temperature  T  =  298°K  ( (3  Is  the 
Volume  Percentage  of  Fuel  In  the  Initial  Mixture) 
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However,  the  error  due  to  this 
fact  Is  small,  since  the  ratio  of 
the  vibrational  and  rotational 
partition  functions  of  the  mole¬ 
cules  and  their  Ions  In  the  range 
of  our  temperatures  ( 1500-4000°K) 
Is  close  to  unity  [8].  Nonethe¬ 
less,  for  the  Ionization  process 
NO  ^  N0+  +  e,  which  makes  the 
principal  contribution  to  the 
electron  concentration  in  the  com¬ 
bustion  products  of  pure  mixtures 
(without  additives),  and  also  for 
the  main  processes  Involving  the 
formation  and  disintegration  of 
negative  .Ions 

0H“  OH  +  e  and  0”  5*  0  +  e 
(see  below)  we  shall  use  hence¬ 
forth  the  available  more  accu¬ 
rately  calculated  equilibrium  con¬ 
stants*  listed  In  Table  3* 

The  electron  concentration 
Is  determined  from  the  formula 

ne  "  Pe1^  U-9) 
Here  pft  and  £  are,  respec¬ 
tively,  the  partial  pressure  of 
the  electrons  and  the  over -all 
pressure  of  the  mixture  [atm]. 


and  n  Is  the  total  concentration 


of  the  molecules  at  a  given  pressure  and 
tempera  tur e  [ cm"^ ] . 

An  analysis  of  Relations  (1.8), 
(1.5)*  and  (1.9)  shows  that  the  concen¬ 
tration  of  the  electrons  of  a  multicom¬ 
ponent  mixture  is  determined  essentially 
by  the  component  having  the  lowest  ion¬ 
ization  potential,  although  it  may  be 
present  in  an  insignificant  amount.  With 
further  increase  in  the  temperature,  when  this  component  turns  out  to 
be  practically  completely  ionized  and  ceases  to  play  a  role  in  the 
total  ionization  of  the  mixture,  the  electron  ionization  is  determined 
successively  by  each  following  more  readily  ionized  remaining  compo¬ 
nent,  up  to  its  complete  ionization.  The  concentration  of  the  negative 
ions,  in  contrast  to  the  ionization  process,  is  determined  essentially 
by  the  component  that  has  the  highest  electron  affinity,  and  with  in¬ 
crease  in  temperature  their  concentration  decreases.  Naturally,  the 
presence  of  electro-negative  components  in  the  mixture,  other  condi¬ 
tions  being  equal,  reduces  the  concentration  of  the  electrons  owing  to 
their  capture  by  these  components. 

It  is  known  that  in  order  to  increase  the  electric  conductivity 
of  the  gases  one  adds  to  them  small  amounts  of  easily  ionized  sub¬ 
stances,  most  frequently  alkali  metals,  which  have  the  lowest  ioniza¬ 
tion  potentials.  It  is  necessary  to  note,  however,  that  when  alkali 
metals  are  added  to  combustible  mixtures  it  is  necessary  in  the  gen¬ 
eral  case  to  take  into  account  the  production  of  compounds  between  the 
alkali  metals  and  certain  components  of  the  combustion  products  and 
the  resultant  concentration  drop  of  these  metals  and  hence  also  the 
reduction  in  the  electron  concentration  in  the  mixture.  The  most  prob- 
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.TABLE  3 

Equilibrium  Constant  for 
the  Processes 


NO  C  NO*  Oil"  COII  +»,0-50  +  « 


T*K 

*NO 

*OH- 

2:v)0 

0.125* 

0.1521 

24oo 

1  .3GG*  10"w 

0.2017 

0.2304 

25**0 

9.012* 10-» 

0.3141 

0.3-189 

2f>00 

5.i62*10~u 

0.4746 

0.4859 

27U0 

_ 

0.0981 

0.6806 

2800 

— 

1 .00* 

0.9V57 

21**0 

— 

1.409 

1.257 

*nno 

— 

1 .94* 

1.665 

*ino 

_ 

2.0*2 

2.170 

*200 

— 

3.505 

2.790 

able  are  the  hydroxides  of  the  alkali  metals  CsOH  [9,  10],  KOH  [11,  12] 
NaOH  [9]  and  the  lower  oxide  of  potassium  KgO  [11];  the  most  stable  of 
these  at  higher  temperatures  are  CsOH  and  KOH. 

Let  us  estimate  approximately  the  amount  of  alkali  metal  A  com¬ 
bining  with  the  hydroxyl.  We  use  for  the  process  AOH  £  A  +  OH  the  re¬ 
lations 


A'aoii  —  ,  Pk  =*  Pa+  Paou 

Paou 


and  assume  that  at  small  amounts  of  the  additives  pQjj  «  PoH*  then 
obtain  the  ratio  of  the  partial  pressure  of  the  bound  metal  (hydroxide) 
pA0H  t0  the  partial  pressure  of  the  metal  vapor  p^ 


AaOH  Poh 

*0H  +  *AOH 


(1.10) 


According  to  [9-10]  for  the  process  CsOH  ;*  Cs  +  OH  and  KOH  K  + 

+  OH  at  a  temperature  T  =  2000°K  the  equilibrium  constants  are,  respec- 

..Jl  O 

tlvely,  Kqb0H  =  3*2*10  and  Kjqjj  =  10  J  atm,  while  the  thermal  ef¬ 
fects  are  91  and  86  kcal/mole.  Extrapolating  these  constants  to  the 
values  of  the  temperatures  of  the  combustible  mixtures  which  we  are 
calculating  (Table  2),  we  obtain  from  Formula  (1.10)  that  In  the  com¬ 
bustion  products  of  fuel-oxygen  mixtures,  owing  to  their  high  tempera¬ 
tures,  there  is  practically  no  cesium  or  potassium  combined  with  the 
hydroxyl,  while  In  the  combustion  products  of  fuel-air  mixtures,  which 
have  a  lower  temperature,  not  more  than  30$  cesium  and  10$  potassium 
combined  with  the  hydroxyl  can  be  found;  to  simplify  the  calculations 
of  the  electron  concentration  and  of  the  electric  conductivity,  this 
will  be  disregarded  from  now  on. 

2.  Initial  data  for  the  calculation  of  the  electron  concentration 
and  the  electric  conductivity  of  gas  mixtures.  The  values  of  the  Ion¬ 
ization  potentials  and  the  ratios  of  the  statistical  weights  of  the 
ground  state,  necessary  for  the  calculation  of  the  components  of  dlf- 


TABLE  .4 

Values  of  the  Electron  Concentration,  n_  in  cm  J, 

and  the  Electric  Conductivity  a  in  ohm"*  cm”  ,  for 
Inert  Gases  with  1 %  by  Weight  of  Cesium,  Potassium, 
and  Sodium  Additive,  and. also  for  Pure  Vapors  of 
Alkali  Metals  at  Different  Temperatures  T®K  and  at 
a  Pressure  p  =»  1  atm 


1 

t*k 

JL 

Cnee* 

1500 

2000 

2500 

2000 

IV* 

46* 

Ar+l%Cs?* 

4.15*10" 

1.93- 10“ 

1.94*10" 

8.65*10" 

2.25*10“ 

3.74*10“ 

0 

5.91  •!<>-» 

2.89*10-* 

2.18 

4.83 

7.70 

10.20 

Ar+l**"* 

1.36*10" 

9.68-10*' 

1.27*10" 

7.04-10** 

2.32*10“ 

5.30*10“ 

a 

6.36-10'* 

5.15*10-* 

6.50-10'* 

2.63 

5.90 

9.43 

N,+l%Cs"* 

3.45*10" 

1.60*10“ 

1.6010“ 

7.90-10*' 

1.80*10“ 

2.83*10" 

a 

1. 19-10'* 

5.57*10'* 

5.06- 10-» 

1.87 

3.86 

5.57 

He 

1. 11-10" 

8.07*10" 

1.06*10" 

5.84*10" 

1.92-10** 

4.2910" 

N,+i%K 

a 

2.86*10”* 

2.19*10”* 

2.7910-* 

1.37 

3.55 

6.30 

He 

6.71*10* 

1.04- 10** 

2.1610" 

1.65*10" 

7.01*10" 

2.05*10" 

N,+i%N» 

1.58-10”* 

2. 61-10'* 

5.54-10-' 

4.19*10”* 

1.59 

3.85 

*  He+l%Cs* 

1.32*10" 

6.07*10" 

5.98-10“ 

2.40*10" 

4.60*10“ 

5.20*10" 

0 

3.84-10”* 

2.02-10”* 

2.12-10-* 

8.40*10-* 

1.59 

1.93 

Its 

4.32*  1010 

3.08*10" 

3.98-10“ 

2.17*10" 

6.41*10" 

9.50*10" 

He+1%* 

a 

1.20-10”4 

9.80* 10-* 

1.38-10*1 

7.45*10-* 

1.97 

2.89 

!00%Cs 

a 

7.58-1011 

3.54*10" 

3.55-10“ 

1.68*10" 

4.97*10" 

1.11*10" 

3.80*10”* 

2.01-10'* 

2.14-10-' 

1.08 

3.09 

6.S7 

n« 

100%  K 

e 

1. 35*10" 

9.59*10" 

1.26-10“ 

7.06*10“ 

2.40-10"] 

6.03*10" 

6.80*10-* 

5.54*10"* 

8.05-10-' 

4.78*10”* 

1.62 

3.99 

100%N«  "* 
a 

6.12*10" 

9.49*10" 

1.97-10“ 

1.51*10“ 

6.46*10“ 

1 

1  jO-10" 

4.1010'* 

7.34*10-* 

1.70-10-' 

1. 40*10-* 

6.18*10-* 

1.81 

l)  Mixture. 


ferent  mixtures  and  products  of  combustion,  are  chosen  in  accord  with 
Table  1,  from  which  it  is  seen  that  the  smallest  ionization  potentials 
are  possessed  by  the  alkali  metals  Cs,  K,  Na,  and  nitric  oxide  NO.* 

As  far  as  we  know,  of  all  the  components  indicated  in  the  table  the 
largest  electron  affinity  is  possessed  by  OH  —  1.735  ev  [22],  0  —  1.465 
ev  [23,  24],  and  Og  -  0.87  ev  [25]. 

The  equilibrium  chemical  composition  of  the  products  and  the  com- 
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bustion  temperature  T#  in  °K  of  the  combustible  mixtures  employed*  are 
listed  in  Table  2. 


The  electric  conductivity  of  the  mixture  is  determined  from  the 


Pig.  1.  Dependence  of  the 
electric  conductivity  on 
the  temperature  for  inert 
gases  (argon,  helium, 
nitrogen)  with  additives 
of  alkali  metals  (cesium, 
potassium,  sodium)  and 
pure  vapors  of  alkali 
metals:  l)  Ar  +  156  Cs;  2) 
Ar  +  1J6  K;  3)  Ng  +  156  Cs; 

4)  Ng  +  lj6  K;  5)  Ng  +  1$6 

Na;  6)  He  +  1J6  Cs;  7) 

He  +  156  K;  8)  lOOjC  Cs;  9) 
IOO56  K;  10)  IOO56  Na. 


relation  [27*  28) 


0.532 


"A 

< 


(2.1) 


where  e  and  m.  are,  respectively,  the 
charge  and  mass  of  the  electron,  n  the 

W 

electron  concentration  [cm-^],  k  Boltz¬ 
mann's  constant,  T  the  temperature  [°K], 
and  n^  and  the  concentration  [cm“^] 
and  the  average  electronic  diffusion  col¬ 
lision  cross  section  for  neutral  par- 

A 

tides  or  ions  of  the  mixture  [cm  ].  For 
ions,  the  latter  is  calculated  from  the 
formula  [28] 


Qh  =  0.C0 


Hiere  are  few  reliable  data  on  the 


cross  sections  far  the  collisions  of  electrons  with  the  neutral  compo¬ 
nents,  pertaining  to  the  slow  electrons  with  energy  less  than  0. 5  ev 
which  are  of  Interest  to  us.  In  the  calculation  of  the  electric  con¬ 


ductivity  of  mixtures  of  argon,  helium,  and  nitrogen  with  additives  of 
cesium,  potassium,  and  sodium  in  the  temperature  interval  from  1500  to 
4000°K,  we  assumed  the  following  values  of  the  cross  sections:  for  Ar  — 
2-10"^  [29,  30],  He  —  6.8*10"^  [31  ],  Ng  —  depending  on  the  tempera¬ 
ture  (5.0-7.4)-10“16  [29],  Cs  and  K-  4*10"14  [32],  NS  -  3*10^  cm2 
[32]. 

Hie  cross  sections  for  the  collisions  between  the  electrons  and 
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the  combust ion -product  components  were  taken  to  be:  for  H  —  1.5*10 
[33],  Hg  —  1.2- 10-15  [32,  34],  0-  (1.6-2.1).10-15,  o2  -  (3-3-3. 6)  x 
x  10-16,  Ng  —  (6.3-6.4)-10'16,  NO  -  (5.4-5.6)*10"16  [29],  CO  -  (7-8)  x 
x  10-16,  C02-  (1.2-3)-10"15  [33],  HgO  —  (5-7)*10-15  [33,  34],  OH - 
2*10-*-*  cm^  (tentative). 

It  follows  from  the  foregoing  data  that  the  largest  collision 
cross  sections  are  possessed  by  the  alkali  metals,  and  the  smallest  by 
argon  (the  Ramsauer -Townsend  effect). 

3,  Results  of  the  calculations,  a.  Electron  concentration  and 
electric  conductivity  of  inert  gases  with  alkali-metal  additives.  The 
calculated  dependence  of  the  electron  concentration  and  of  the  elec¬ 
tric  conductivity  of  the  mixture  on  the  temperature,  pressure,  and 
amount  of  additive  has  been  determined  in  the  range  of  variation  of 
practical  interest  for  each  of  these  parameters. 

The  results  of  calculation,  based  on  Relations  (1.5),  (1.9),  and 
(2.1),  of  the  electron  concentration  and  of  the  electric  conductivity 
of  the  inert  gases  argon,  helium,  and  nitrogen*  with  additives  of  cesium 
and  potassium,  and  in  the  case  of  nitrogen  also  with  additive  of  sodium, 
in  the  temperature  interval  from  1500  to  4000°K  at  atmospheric  pres¬ 
sure  and  at  a  constant  additive  content  amounting  to  1 %  of  the  weight 
of  the  mixture  are  listed  in  Table  4,  while  the  data  for  the  electric 
conductivity  are  plotted  in  Pigs.  1  and  2.  For  the  sake  of  comparison 
we  show  here  also  the  data  calculated  under  these  conditions  from  For¬ 
mulas  (1.4),  (1.9),  and  (2.1)  for  the  electric  conductivity  of  the 
pure  cesium,  potassium,  and  sodium  vapor. 

It  follows  from  the  curves  of  Fig.  1  that  in  the  entire  tempera¬ 
ture  interval,  the  largest  electric  conductivity,  reaching  a  value  a  « 

~  10  ohm-'*'  cm-*  at  a  temperature  T  =  4000°K,  is  possessed  by  argon 
with  cesium  additive.  In  the  central  part  of  the  calculated  tempera - 
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Fig.  2.  Dependence  of  the  elec¬ 
tric  conductivity  on  the  recip¬ 
rocal  temperature  for  inert 
gases  (argon,  helium,  nitrogen) 
with  additives  of  alkali  metals 
(cesium,  potassium,  sodium)  and 
pure  vapors  of  alkali  metals: 

1)  Ar  +  1%  Cs;  2)  Ar  +  1%  K; 

3)  N2  +  1 %  Cs;  4)  N2  +  1*  K; 

5)  N2  +  1%  Na;  6)  He  +  1$6  Cs; 

7)  He  +  ljf  K;  8)  IOC#  Cs;  9) 

IOC#  K;  10)  10C#  NS. 

ture  interval,  approximately  between  2000  and  3300°K,  the  electric  con 

ductivlty  of  argon,  nitrogen,  and  helium  is  successively  decreased  rel 

atlve  to  the  value  with  cesium  and  potassium  additives.  The  smallest 

electric  conductivity  is  in  this  case  that  of  the  mixture  of  nitrogen 

with  sodium. 

The  electric  conductivity  of  IOC#  sodium,  potassium,  and  cesium 
vapor  (see  the  dashed  lines  on  the  figures)  rise  with  decreasing  ion¬ 
ization  potential  from  sodium  to  cesium,  which  is  natural,  but  with 
the  exception  of  a  small  region  of  low  and  high  temperatures,  the  elec 
trie  conductivity  of  10C#  cesium  vapor  turns  out  to  be  lower  than  the 
electric  conductivity  of  argon  and  nitrogen  with  cesium  and  potassium 
additives,  while  the  electric  conductivity  of  10C#  potassium  vapor  is 
lower  than  the  electric  conductivity  of  helium  mixed  with  potassium 
and  cesium,  owing  to  the  small  mean  free  path  of  the  electrons  due 


to  the  larger  cross  sections  for  the  collision  between  the  latter  and 
the  alkali  metalB. 

As  can  be  seen  from  Fig.  2,  the  temperature  variation  of  the  elec¬ 
tric  conductivity  at  low  degree  of  ionization  of  the  mixture  (up  to  a 
temperature  ~2500°K  for  the  argon-cesium  mixture  and  up  to  a  tempera¬ 
ture  «3000°K  for  the  remaining  mixtures)  can  be  approximated  with  suf¬ 
ficient  accuracy  by  an  exponential  dependence  of  the  type  a  = 

=  C  exp  (— D/T).  With  further  increase  in  the  temperature,  and  conse¬ 
quently  also  in  the  ionization,  the  increase  in  the  electric  conduc- 
t_vity  slows  down  more  and  more  as  the  initial  substances  become  con¬ 
sumed  and  because  of  the  relative  decrease  in  the  mean  free  path  of 
the  electrons,  owing  to  the  more  probable  collision  between  them  and 
the  ions  (although  the  electron-ion  collision  cross  section  decreases 
with  increasing  temperature,  nevertheless  it  is  several  orders  of  mag¬ 
nitude  larger  than  the  cross  section  for  neutral  molecules). 

The  calculated  data,  in  the  same  temperature  interval,  for  the 
dependence  of  the  electron  concentration  and  the  electric  conductivity 
of  the  mixture  on  the  pressure,  in  its  range  of  variation  from  0.01  to 
10  atm,  using  a  mixture  of  nitrogen  with  1; %  potassium  by  weight  as  an 
example,  are  listed  in  Table  5  and  are  plotted  for  the  electric  con¬ 
ductivity  in  Fig.  3-  It  follows  from  these  data  that  the  electric  con¬ 
ductivity  of  the  mixture  decreases  with  increasing  pressure,  and  up  to 
a  temperature  s=2500°K  the  Increase  follows  the  law  =  const,  i.e. , 

it  is  Inversely  proportional  to  the  square  root  of  the  pressure. 

In  the  general  case,  for  a  small  degree  of  ionization  (on  the  or¬ 
der  of  3$6),  we  can  assume  approximately  that  the  dependence  of  the 
electric  conductivity  of  the  given  mixture  on  the  pressure  and  on  the 
temperature  satisfies  the  equation 
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in  particular,  for  a  mixture  of  nitrogen  with  1%  potassium  by  weight 
we  have 

.  •tt-W* /  —2.6- 10*  \ 

* "  \  f — ) 


The  dependence  of  the  electric  conductivity  of  the  above-mentioned 
mixtures  of  inert  gases  with  additives  of  alkali  metals  on  their  weight 
fraction  y  in  the  mixture,  at  atmospheric  pressure  and  at  a  temperature 
T  =  3000°K,  is  shown  in  Pig.  4. 


TABLE  5 

Values  of  the  Electron  Concentra¬ 
tion  n_,  in  cm"^,  and  of  the 
Electric  Conductivity  o,  111 

ohm”1  cm”1,  for  Nitrogen  with  1% 
Potassium  by  Weight  Added  at  Dif¬ 
ferent  Pressures  £  in  atm  and 
Temperatures  T  in  °K 


rx 

P* 

Am.1 

Ml 

•.1 

1  “ 

i  • 

1500 

I.I1-10** 

3.57.10* 

1.1110” 

3.57*19* 

2000 

8.07- 1011 

2.82.10“ 

8.0710“ 

2.36-M* 

2500 

1 .03.10“ 

3.35.10* 

1 .06-10“ 

3.3510* 

:mm)o 

4.94.10** 

1.78-10* 

5.84-10* 

1.88-49* 

5500 

I.UMO* 

5.21-10“ 

1.9210“ 

8.36-10* 

4000 

1.30-10* 

9.6510* 

4.29-10“ 

t. 56-49* 

1500 

2.80-10* 

9.22-10* 

2. 86-10"* 

9.24-40* 

o 

2000 

2.01-10* 

6.80.10 r* 

2.19-10* 

7.00-40* 

2500 

1.53 

7.5210* 

2.79-10* 

9:49-40* 

5000 

3.69 

2.70 

1.37 

5.49-10* 

3500 

5.30 

5.20 

3.55 

4.87 

4000 

6.48 

7.00 

6.30 

3.75 

1)  Atm. 


Prom  the  curves  in  this  figure,  which  were  calculated  from  For¬ 
mula  (2.1),  we  see  that  for  specified  conditions  with  respect  to  pres¬ 
sure  and  temperature,  a  maximum  electric  conductivity  of  the  gas  Is 
attained  at  a  certain  (optimal)  amount  of  additive.  Further  increase 
of  the  additive  leads  to  a  reduction  in  the  electric  conductivity. 
Qualitatively  this  follows  directly  from  Formula  (2.1),  according  to 
which  the  electric  conductivity  increases  with  increasing  electron  con- 
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Pig.  3  Pig.  4 


Fig.  3.  Dependence  of  the  electric  conductivity  of 
nitrogen  +  1%  potassium  b’r  weight  on  the  temperature 
at  different  pressures:  1)  0.01  atm;  2)  0.1  atm;  3) 

1  atm;  4)  10  atm. 

Fig.  4.  Dependence  of  the  electric  conductivity  of 
inert  gases  (argon,  helium,  nitrogen)  with  alkali 
metals  added  (cesium,  potassium,  sodium)  on  the 
weight  fraction  of  the  additives  (T  =  3000°K,  p  =  1 
atm):  l)  Ar  +  Cs;  2)  Ar  +  K;  3)  Ng  +  Cs;  4)  Ng  +  K; 

5)  Ng  +  Na;  6)  He  +  Cs;  7)  He  +  K. 

centration  and  decreases  with  increasing  electronic  collision  cross 
section.  Naturally,  the  increase  in  the  amount  of  Impurity  in  the  form 
of  alkali  metals,  which  have  much  higher  cross  sections  compared  with 
ordinary  gases,  leads  on  the  one  hand  to  an  increase  in  the  electric 
conductivity  due  to  the  increase  in  the  electron  concentration  in  the 
mixture,  and  on  the  other  hand  to  a  decrease  due  to  the  increase  in 
the  total  electronic  collision  cross  section: 

0  - 

i 

As  a  net  result,  at  a  definite  (optimal)  amount  of  impurity,  a 
maximum  of  electric  conductivity  sets  in,  and  the  dependence  of  the 

maximum  on  various  factors  can  be  established  by  an  analysis  of  the 

indicated  formula. 

For  the  sake  of  simplicity  let  us  take  a  two-component  mixture, 
consisting  of  a  fundamental  component  of  some  electro-positive  gas  and 
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the  vapor  of  some  additive,  and  neglecting  the  ionization  of  the  .fun¬ 
damental  component  let  us  assume  that  In  our  temperature  range  only 
the  additive  substance  Is  ionized.  Expressing  in  Formula  (2.1)  the  ad¬ 
ditive  concentration  n^,  the  concentration  of  the  main  component  ng, 
and  the  electron  concentration  n  in  terms  of  the  corresponding  par- 
tial  pressures  p1#  pg,  and  p@,  we  obtain  with  the  aid  of  (1.1),  (1.3), 
and  (1.9) 

•  _  0 .3.W 

Differentiating  this  expression  with  respect  to  p_  under  the  aim- 
pllfylng  assumption  that  the  collision  cross  section  of  the  electrons 
with  the  positive  ions  of  the  additive  q£  does  not  depend  on  the  par¬ 
tial  pressure  pe  of  the  electrons  (the  error  introduced  by  this  in  the 
electric  conductivity  is  insignificant  in  our  case),  and  elating  the 
first  derivative  to  zero,  we  obtain  a  relation  for  the  optimum  partial 
pressure  of  the  electrons: 

Using  this  equation  along  with  Eqs.  ( 1. 1) -( 1. 3)  for  the  ioniza¬ 
tion  of  the  additive,  we  obtain  the  optimum  ratio  aQpt  of  the  weight 
of  the  additive  0^  to  the  weight  of  the  gas  Og 


where  and  |Xg  are,  respectively,  the  molecular  weight  of  the  addi¬ 
tive  and  the  gas. 

If  we  assume  that  the  collision  cross  sections  are  Independent  of 
the  temperature,  then  it  follows  from  this  formula  that  the  weight  ra¬ 
tio  Ogpt  increases  with  increasing  temperature,  owing  to  the  increase 
in  the  ionization  constant  K^,  and  decreases  with  the  increase  in  the 
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total  pressure  £, 

Neglecting  in  the  numerator  of  (3.1)  the  expression  under  the 
square  root  and  neglecting  In  the  denominator  both  terms  compared  with 

unity,  we  can  assume  for  a  clearly  ten¬ 
tative  estimate  under  the  conditions  con¬ 
sidered 

„  _  i*i 

a°*  “ 

which  coincides  in  the  particular  case 
when  Qg  «  with  the  concentration  ra¬ 
tio  [351. 

The  values  of  the  optimum  weight 
fraction  of  the  additive  in  the  mixture 

~  .  —  “"P* 

Top*  -  i  +  «opt 

and  of  the  corresponding  maximum  electric  conductivity  °max  at  atmos¬ 
pheric  pressure  and  at  a  temperature  T  =  3000°K  are  listed  in  Table  6 
for  various  mixtures.  It  is  seen  from  this  table  that  for'  a  mixture  of 
argon  with  potassium  and  cesium,  only  a  very  insignificant  amount  of 
optimal  additive  (less  than  1#)  is  necessary  to  attain  maximum  elec¬ 
tric  conductivity  of  3*57  and  5-2  ohm”1  cm“\  respectively;  for  mix¬ 
tures  of  nitrogen  with  sodium,  potassium,  and  cesium  this  value  does 
not  exceed  8#,  and  for  mixtures  of  helium  with  potassium  and  cesium  it 
reaches  37#  at  relatively  lower  absolute  values  of  the  maximum  elec¬ 
tric  conductivity. 

b.  Concentration  of  electrons  and  electric  conductivity  of  com¬ 
bustion  products.  The  concentration  of  the  electrons,  determined  on 
the  basis  of  the  Saha  equation  (1.8),  and  the  electric  conductivity  of 
combustion  products  without  additives  of  readily  ionized  substances 
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ni/n2  = 


TABLE  6  . 

Values  of  the  Optimum 
Weight  Fraction  7opt  # 

and  the  Maximum  Electric 
Conductivity  omax, 

ohm-1  cm”1,  for  Different 
Mixtures  at  Pressure  p  = 
=  1  atm  and  Temperature 
T  =  3000°K 


1  CwCfc  |  Y0pt  | 

•msi 


Ar  -f  Cs 

0.22 

5.20 

Ar-f  K 

O.UO 

3.57 

Ho+G» 

37. 00 

2.71 

He-f  K 

14.80 

1.50 

N,  +  Oi 

7.SJ5 

2.71 

N,+  K 

2.42 

1.50 

1.85 

0.44 

l)  Mixture. 


are  small  and  according  to  calculation,  with  allowance  for  the  data  of 
Table  3,  do  not  exceed*  ng  =  1.86*1010  cm”3  and  o  =  3.53-10”5  ohm”1  cm"1, 
the  latter  figures  pertaining  to  the  combustion  products  of  acetylene 
with  03Qrgen.  Addition  of  easily  ionized  substances,  particularly  cesium 
and  potassium,  in  amounts  of  Vf>  by  weight,  as  can  be  seen  from  Table  7# 
greatly  increases  the  electron  concentration  and  the  electric  conduc¬ 
tivity  of  the  combustion  products. 

TABUS  7 

Values  of  the  Electron  Concentration  n„,  in  cm”3 

e  -1  -l 

and  of  the  Electric  Conductivity  a  in  ohm  cm 
for.  Combustion  Products  with  1%  by  Weight  of 
Cesium,  Potassium,  and  Sodium  Added  at  a  Pressure 
pal  atm 


1  r»pmm  mmok 

c 

% 

K 

N. 

% 

0. 

*• 

** 

• 

2  Boxopoa  — 

HBRH 

0.47-lfr* 

O  BOUTS 

0  IWwOfWMI  — 

3005 

5.62. 10“ 

3.64*  10*1 

4.92-10“ 

3.1310-1 

1.45-10* 

9.44-10“* 

il  KBOJIOfMJI 

“  Awnuwfi  — 

2575 

1.98*10“ 

3.52-10"* 

1. 36.10* 

2.19*10-* 

2.93*10“ 

4.69- 10"* 

c  .  swAir* 

J  A  firm  mo  — 

3210 

1.04-10“ 

1.50 

9.55-10“ 

1.27 

3.02*10“ 

4-34-10"* 

mwflopoji 

O  Okucii  yrMpou  — 
0*0*0 pox 

2965 

7.18-10“ 

1.20 

5.81-10“ 

1.02 

1.62*10“ 

2.97-10"* 

1)  Combustible  mixture;  2)  hydrogen-air;  3)  hydro- 
gen-oxygen;  4)  acetylene-air;  3)  acetylene -oxygen; 
6)  carbon  monoxide -oxygen. 


The  greatest  value  of  the  electron  concentration  and  the  electric 
conductivity  is  possessed  in  this  case  (Table  7)  by  the  combustion 
products  of  acetylene  and  of  carbon  monoxide  with  oxygen.  It  must  be 
noted  that  in  spite  of  the  fact  that  the  temperature  of  the  combustion 
products  of  hydrogen  with  oxygen  is  100°  higher  than  the  tempera ture 
for  carbon  monoxide  with  oxygen,  however,  for  an  almost  equal  electron 
concentration,  their  electric  conductivity  (owing  to  the  larger  value 
of  the  collision  cross  section  Q»  of  the  combustion  products  (Table  8) 
of  this  mixture  the  main  contribution  is  made  by  the  water)  is  several 
times  smaller  than  the  electric  conductivity  of  the  combustion  products 


TABLE.  8 

Values  of  Optimal  Content  V0pt  of  Additives  of 

Cesium,  Potassium,  and  Sodium  and  of  the  Maximum 

Electric  Conductivity  omnT  In  ohm”*  cm"*  for  the 

Combustion  Products  of  Combustible  Mixtures  with 
Alkali-Metal  Additives 


1  ropwwian 
cnee* 

CJ*» 

c« 

K 

VK 

N 

T©i* 

•max  1 

"mat 

*©P4 

"mat 

pBojtnpfti  — 

2387 

24.3 

2.7 

30.0 

0.27 

11. 0 

0.10 

9.5 

0.0 

0  *MAyt 

0.79 

3065 

14.6 

4.2 

54.0 

1.83 

26.7 

.  24.0 

0.29 

1  HUC 

l}AuerM/ieN  — 

^  Ann  mew  — 

2575 

25.2 

1.3 

15.6 

0.74 

5.2 

0.30 

4.2 

0.06 

32t0 

26.3 

ft  .5 

18.4 

3.38 

6.0 

1.85 

4.9 

0.60 

KMOAOpoa 

i  Owwck  yr.iepo 

«Ut  -  MNCXOpQJ 

2.9 

2965 

35.4 

1.0 

9.3 

2.36 

1.41 

2.3 

0.33 

l)  Combustible  mixture;  2)  hydrogen-air;  3)  hydro¬ 
gen-oxygen;  4)  acetylene -air;  5)  acetylene -oxygen; 
6)  carbon  monoxide -oxygen. 


of  carbon  monoxide  with  oxygen.  The  electric  conductivity  of  combus¬ 
tion  products  with  potassium  additive  is  close  to  that  with  an  addi¬ 
tive  consisting  of  the  same  fraction  by  weight  of  cesium  (the  differ¬ 
ence  is  approximately  by  1.5  times),  owing  to  the  small  difference  In 
the  ionization  potentials  and  the  relatively  larger  concentration  of 
the  potassium  particles  per  unit  volume  as  compared  with  cesium,  re¬ 
sulting  from  the  much  smaller  molecular  weight  (p.R  =  39.1,  =  132.9). 

Although  the  data  of  Table  7  have  been  obtained  with  account  of  the 
formation  of  negative  ions,  the  latter  can  be  neglected  at  high  tem¬ 
peratures  of  the  combustion  products,  for  even  under  the  most  favor¬ 
able  case,  at  a  maximum  amount  of  the  hydroxyl  OH  and  atomic  oxygen 
(hydrogen -oxygen  and  acetylene-oxygen  combustible  mixture)  this  leads 
only  to  an  error  that  does  not  exceed  5#  in  the  direction  of.  higher 
electron  concentration. 

If  we  regard  the  combustion  products  of  each  of  the  combustible 
mixtures  as  a  homogeneous  gas  having  the  same  molecular  weight  as  n# 
and  total  cross  section  Q*  as  the  combustion  products  of  the  mixture 
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(i.e.,  ii  =  p.*»  .Q  =  Q*)»  then  according  to  calculation  by  means  of  For¬ 
mulas  (3*1)  and  (2.1)  we  obtain  the  values  of  the  optimal  additive  of 
cesium,  potassium,  and  sodium,  t,  and  the  corresponding  maximum 
electric  conductivity  omaT  of  the  combustion  products  of  the  combus¬ 
tible  mixtures  with  additive,  as  listed  In  Table  8.  What  Is  striking 
here  is  the  high  value  of  the  optimal  additive  of  cesium  (up  to  54j6) 
for  the  combustion  products  of  all  the  considered  combustible  mixtures, 
and  also  that  of  potassium  (up  to  26.7#)  and  sodium  (up  to  24 #)  for 
the  combustion  products  of  hydrogen -oxygen  and  hydrogen-air  mixtures^ 
For  the  remaining  combustible  mixtures,  the  value  of  the  additive 
ranges  between  2.3  and  6#,  which  is  probably  tolerable  in  practice, 
although  it  should  reduce  somewhat  their  combustion  temperature.  Hie 
value  of  the  electric  conductivity,  however,  does  not  exceed  in  this 
case  I.85  ohm”^  cm”'*'  (the  combustion  products  of  the  acetylene -oxygen 
mixture  with  6#  potassium),  so  that  the  burned  gases,  even  when  al¬ 
kali  metals  are  added,  are  unfortunately  relatively  poor  electric  con¬ 
ductors. 

In  conclusion  it  must  be  noted  that  the  procedure  of  adding  easily 
ionized  substances  to  the  working  medium  is  essentially  a  method  of 
reducing  the  temperature  of  the  most  difficultly  ionized  working  medium 
when  a  definite  electric  conductivity  is  attained.  Successful  searches 
for  methods  of  increasing  the  electric  conductivity  and  further  de¬ 
creasing  the  temperature  of  the  working  medium  to  sensible  limits  with¬ 
out  reducing  the  electric  conductivity,  along  with  the  development  of 
new  high- temperature  materials,  are  exceedingly  important,  since  they 
would  facilitate  in  principle  the  solution  of  many  applied  problems  of 
practical  Importance,  particularly  problems  connected  with  magnetogas- 
dynamic  conversion  of  heat  into  electricity. 
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[Footnotes  ] 


Even  for  such  easily  ionized  substances  as  alkali  metals, 
the  error  resulting  from  this  is  negligibly  small  and  there 
is  no  sense  in  taking  a  practical  account  of  it. 

The  reduction  in  the  ionization  potential  [5,  7]  resulting 
from  the  screening  of  the  particle  by  the  fields  of  the  sur¬ 
rounding  charges  is  disregarded. 

The  data  were  graciously  furnished  to  us  by  L.  V.  Gurevich. 

The  ratio  of  the  statistical  weights  2<Jq  /q^  in  the  calcu¬ 
lation  of  the  equilibrium  constant  of  the  process  of  forma¬ 
tion  end  disintegration  of  negative  ions  of  molecular  oxygen 

is  assumed  to  be  1.5. 

The  data  for  the  hydrogen-oxygen  mixture  were  borrowed  from 
[26],  and  the  data  for  the  remaining  combustible  mixtures 
were  graciously  furnished  by  I.V.  Veyts. 

We  disregard  dissociation  of  molecular  nitrogen. 

We  note  that  the  actually  measured  concentration  of  the  elec¬ 
trons  in  the  flame  zone  itself  is  two  to  four  orders  of  mag¬ 
nitude  larger,  depending  on  the  composition  of  the  combus¬ 
tible  mixture,  than  that  thermodynamically  calculated  with 
the  Saha  equation  [3,  22,  36].  It  is  easy  to  see,  however, 
that  in  this  case  the  electric  conductivity  of  the  flame 
still  remains  low. 
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INFLUENCE  OF  THE  DIFFUSION  OF  RADIATING  IMPURITIES 
ON  CONVECTIVE  HEAT  EXCHANGE 

A.N.  Rumynskly 
(Moscow) 

We  consider  the  boundary  layer  near  the  forward  critical  point  of 
an  axially  symmetrical  blunt-nosed  body,  from  the  surface  of  which  ac¬ 
tive  (radiating)  gas  Is  released.  It  Is  assumed  that  the  boundary 
layer  represents  a  binary  mixture.  In  which  no  chemical  reactions  take 
place.  To  separate  the  effect  of  the  influence  of  the  diffusion  of 
the  radiating  gas  on  the  heat  exchange,  the  influence  of  the  radiation 
of  the  gas  components  of  the  main  stream  on  the  flow  of  the  gas  is 
neglected.  Under  the  indicated  assumptions,  and  using  a  few  supplemen¬ 
tary  simplifying  premises,  formulas  are  derived  for  the  calculation  of 
the  convective  heat  flows  In  the  vicinity  of  the  critical  point  of  a 
blunted  body  with  radiating  gas  released  from  Its  surface. 

*  *  • 

1.  In  solving  the  problem  we  shall  start  with  the  results  of 
Reference  [l].  In  which  is  considered  heat  exchange  In  a  frontal  point 
over  which  a  radiating  medium  flows,  without  any  gas  released  from  the 
surface  of  the  body.  Since  we  neglect  here  the  influence  of  the  radia¬ 
tion  of  the  main  stream  on  the  flow  of  the  gas  In  the  perturbed  region, 
unlike  the  cited  paper,  there  is  no  need  in  the  instance  under  consid¬ 
eration  to  use  the  complete  system  of  Navler -Stokes  equations,  and  it 
is  sufficient  to  confine  onself  to  an  analysis  of  the  flow  in  a  vis¬ 
cous  boundary  layer  only.  We  neglect  here  the  interaction  of  the  vor- 
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tices  of  the  shock  and  boundary  layers,  i.e.,  we  shall  consider  the 
parameters  on  the  outer  boundary  of  the  boundary  layer  to  be  equal  to 
the  parameters  behind  the  departed  shock  wave,  which  enables  us  to  use 
the  scheme  of  asymptotic  boundary  layer  [2,  3].  The  latter  will  obvi¬ 
ously  be  valid  in  those  cases  when  the  thickness  of  the  boundary  layer 
is  considerably  smaller  than  the  distance  6  between  the  body  and  the 
shock  wave  [4],  as  will  be  assumed  from  now  on. 

The  equations  of  motion  and  of  continuity  for  a  mixture  of  gases 
in  a  boundary  layer,  in  the  presence  of  diffusion  of  the  radiating  im¬ 
purities  released  from  the  surface  of  the  body,  have  the  same  form  as 
in  [1],  and  by  virtue  of  the  statements  made  above  we  must  equate  the 
terms  in  these  equations  containing  the  transverse  pressure  gradient 
to  zero,  and  retain  only  those  terms  which  are  essential  in  the  bound¬ 
ary  layer. 

The  energy  equation  with  account  of  the  heat  flow  transported  by 
diffusion,  heat  conduction,  and  radiation  and  neglecting  the  diffusion 
heat  conductivity,  can  be  represented  in  the  following  form  [5] 

Oi  ifj  d  '  u  9i  \  1 9m  \s  41. 

PU  ~9x  .>jf  oy  \  P  9 y) 

*  P9m  [  P  pP~  PD»  <‘*  ~  **>  -fr]  ( 1- 1 ) 

Here  x  and  %  are  the  distances,  measured  from  the  forward  critical 
point  along  and  normal  to  the  surface  of  the  body,  respectively;  u  and 
v  are  the  velocity  components  along  the  direction;  £  is  the  pressure 
and  q^  the  velocity  of  the  radiant  energy  heat  influx;  the  subscripts 
1  and  2  pertain  to  the  parameters  of  the  active  gas  and  the  main 
stream,  respectively,  while  the  index  s  pertains  to  the  parameters  be¬ 
hind  the  shock  wave;  p,  i,  and  p.  are  the  density,  enthalpy,  and  vis¬ 
cosity  coefficient  of  the  mixture;  P  and  the  molecular  and  diffu¬ 
sion  Prandtl  numbers;  D  is  the  diffusion  coefficient 
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where  ij  and  cpj  are  the  enthalpy  and  the  specific  heat  of  the  J-th 
component,  X  the  heat  conduction  coefficient,  cp  the  average  specific 
heat  of  the  mixture,  T  the  temperature,  lj  the  production  enthalpy  ex¬ 
trapolated  to  absolute  zero,  and  c j  the  concentration  by  weight  of  the 
j-th  component,  satisfying  the  following  continuity  equation: 

P»  JT4  </■=».  2),  Cl=.l-  Ct  (1.2) 

If  we  neglect  thermal  diffusion  for  the  binary  mixture,  the  dif¬ 
fusion  flux  Kj  contained  In  the  right  half  of  (1.2)  is  equal  to 

K,~-Kt~-pDlt£  (1.3) 

In  the  formulation  of  the  problem  as  considered  here,  the  bound¬ 
ary  conditions  for  the  solution  of  the  reduced  system  of  equations- of 
motion,  energy,  and  diffusion  In  the  boundary  layer  will  be 


ib(l,  o  —  r„  i  ■  id  e,  ■  Cj,  for 

■  ct-0  for*--  (1.^) 

The  subscript  w  denotes  here  the  parameters  on  the  wall. 

The  general  expression  for  the  rate  of  radiant- energy  heat  Influx 
Is  given  In  [1,  6]  and  Is  not  presented  here  In  order  to  save  space. 

We  note  only  that  In  the  case  considered  here  the  absorption  coeffi¬ 
cients  of  the  mixture  depend  on  the  concentration  c^  of  the  active  gas, 
with  av  -*•  0  as  c^  -*•  0,  i.e. ,  by  virtue  of  (1.4)  we  have  -*•  0  as  y  -*■  •. 

2.  Using  the  customary  assumption  made  for  the  vicinity  of  the 
critical  point,  namely  that  terms  of  order  0(x  )  are  email,  we  seek 
the  solution  of  the  obtained  system  In  a  form  analogous  to  the  solu¬ 
tion  of  [1) 


«=/’(»!)*.  * 

.  £  =  -  P  v. 


=  — VW— .  i-i(n). 


ci  «  «i  (n) 


(2.1) 
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Here  R  is  the  radius  of  the  midsection  of  the  body,  ub,  are 
the  velocity  and  density  of  the  incoming  stream,  and  b  is  a  constant 
that  depends  on  the  shape  of  the  body  [71. 

The  problem  then  reduces  to  solution  of  the  following  system  of 
ordinary  differential  equations 

v.(/n'  +  2/r-r-^ 

v.  (-y-f) +2/f  +  - - Vft-Wj  “  —  f*  (2.2) 

v.  (£«.’)'  + W--0 


Here 


v. 


1=12-, 

Mb  * 


The  third  term  of  the  left  part  of  the  energy  equation  in  the  sys¬ 
tem  (2.2)  takes  account  of  the  fact  that  when  PD  ^  P  the  energy  trans¬ 
port  processes  occur  at  different  intensities,  owing  to  diffusion  and 
heat  conduction.  When  P^  =  P  this  equation  coincides  in  form  with  the 
equation  considered  in  [1].  The  last  simplification  is  used  to  calcu¬ 
late  the  increment  of  the  convective  heat  flow,  due  to  the  radiation 
of  the  gas  released  from  the  surface  of  the  body.  In  addition,  we  have 
l  =  const.  The  errors  connected  with  the  assumptions  made,  and  the  ad¬ 
visability  of  these  assumptions  in  calculating  the  increment  of  the 
convective  stream,  will  be  discussed  in  the  last  section  of  this  paper. 

Under  the  assumptions  made,  the  system  (2.2)  can  be  written  in 
dimensionless  coordinates  in  the  form 


r+  2 tr  «  r-  j.«'+2/j>r 

C*+2//*,'=  0.  K'=^--±-K,  K  -  (2.3) 

*•  -  *  -  (±JS  fc  -  /.-(P*n)-*/r  9m- 


with  boundary  conditions 
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f  «  0,  /  «  /.,  i  «=  »»,  e,  «=  e,w  for  *  -  • 

/— 1.£— 1,  fl  — 0  for,-..  (2.4) 

The  dimensionless  quantities  are  designated  by  the  subscript  0, 
which  for  simplicity  has  been  left  out  of  the  system  (2.3)-(2.4)j  the 
subscript  m  pertains  to  the  characteristic  mean  values  of  the  corres¬ 
ponding  quantities. 

The  parameter  K*  characterizes  the  main  Influence  of  the  radia¬ 
tion  on  the  characteristics  of  the  boundary  layer  In  the  vicinity  of 
the  critical  point.  Rie  relation  given  above  between  K*  and  the  corres¬ 
ponding  parameter  K,  used  In  [6]  for  a  flat  plate,  characterizes  the 
influence  of  the  pressure  gradient.  The  first  factor  at  K  is  equal  in 
practice  to  1/2  and  appears  only  because  of  some  difference  between 
the  dimensionless  quantities  employed  here  and  the  corresponding  quan¬ 
tities  of  [6].  The  second  factor  is  due  to  the  pressure  gradient  in 
the  vicinity  of  the  critical  point.  With  Increasing  pressure  gradient, 
the  velocity  increases  and  the  time  that  the  particle  stays  in  the 
vicinity  of  the  retardation  point  with  high  values  of  pressure  and  tem¬ 
perature  decreases,  so  that  the  influence  of  the  radiation  diminishes, 
l.e.,  the  parameter  K*  decreases. 

We  shall  assume  that  q^  Is  a  known  quantity,  determined  from  the 
profiles  of  the  temperature  and  of  the  concentration  of  the  active  gas 
in  the  nonradiating  boundary  layer,  the  latter  of  which  is  determined 
by  the  solution  of  the  third  equation  of  the  system  (2.2) 

c* *=  «*.  (*  —  .  AM  =$«*p(—  (2.5) 

•  • 

Here  and  throughout  the  function  f  is  determined  by  the  solution 
of  the  system  (2.2)-(2.4)  with  K*  *  0. 

Taking  Into  account  the  radiation  of  the  active  gas  released  from 
the  surface  of  the  body,  the  distribution  of  the  enthalpy  of  the  bound- 
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ary  layer  Is  written  In  the  form 

<(«-<•+ CAM 

/»(n)  -  jc*p(-2/,j/rfn)j«(xp(2/>j/^)ffa(r,eI)rfnrfn 


(2.6) 


The  heat  flux  transported  to  the  surface  of  the  body  by  diffusion 
and  heat  conduction  Is  determined  by  the  gradients  of  the  temperature 
and  concentration  at  tj  =  0  In  accordance  with  the  formula 

—  »4ir)(,|“**)  (2.7) 


Calculating  the  value  of  the  derivative  of  the  temperature  on  the 
wall  and  substituting  It  into  the  expression  for  the  heat  flux  (2.7) 
we  obtain 


-.)  (2. 8) 

It  was  pointed  out  in  [1]  that  the  variability  of  the  term  con¬ 
taining  the  pressure  gradient  in  the  first  equation  of  the  system  (2.3) 
influences  the  solution  of  the  equation  slightly  and  does  not  affect 
in  practice  the  solution  of  the  energy  equation,  so  that  these  equa¬ 
tions  can  be  solved  independently  of  each  other.  Moreover,  Lees  has 
shown  that  the  right  half  of  the  equation  of  motion  can  be  left  out 
and  only  small  errors  result  from  this,  particularly  at  low  enthalpy 
factors  1^,/lg  «  1,  which  is  usually  the  situation  in  the  cases  of 
greatest  practical  interest.  This  enables  us  to  employ  in  the  numerical 
calculations  the  solution  of  the  Blaslus  equation  for  fw  0,  which 
has  been  tabulated  for  many  values  of  f  in  [8];  in  this  case 


A (1)  -  /, (H)  =  [rp[r-'i* (T,  c,)*)*!  ( 2. 9) 

•  Si 

3.  To  obtain  the  results  in  final  form  and  for  a  quantitative  es¬ 
timate  of  the  effect,  we  shall  assume  the  active  gas  to  be  a  gray  ra¬ 
diator  with  a  volume  absorption  coefficient  k  =  ap,  which  depends  lln- 
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early  on  the  partial  pressure  of  the  active  gas,*  and  we  shall  also 
neglect  the  reabsorption  In  the  boundary  layer.  Under  these  assump¬ 
tions,  the  influx  rate  of  radiant  energy  can  be  represented  In  accord¬ 
ance  with  [1,  6]  in  the  form 

J*  =  Jr*  y(*»T»+  (2  —  *»)  7r*)J  7r“7r»®  (3*1) 

Here  is  the  blackness  coefficient  of  the  surface  about  which 
the  flow  takes  place,  q^  is  the  radiant  flux  incident  on  the  outer 
limit  of  the  boundary  layer  from  the  high-temperature  shock  layer,  and 
o  is  the  Stefan-Boltzmann  constant. 

Since  6/R  «  1  for  hypersonic  velocities,  the  portion  of  the 
shock  wave  that  radiates  in  the  vicinity  of -the  critical  point  (spher¬ 
ical  segment)  can  be  regarded  as  a  plane -parallel  layer  with  thickness 
6,  and  q^,  can  be  determined  from  the  formula 

m  m 

Here  By  is  the  Planck  function. 

If  we  calculate  (jj,  by  using  the  approximate  egression  for  the 
function  E^t)  given  in  (3*2),  then  the  computation  formula  assumes 
the  form 

’  m  •  • 

1r-*'(P.,T:L,)oTt*  A.-P*)  (3*3) 

Here  ee  is  the  effective  blackness  coefficient  of  a  hemispherical 
gas  volume,  radiating  in  the  center  of  the  base  of  the  hemisphere;  I»e 
is  the  equivalent  radius  of  the  considered  radiating  volume.  At  an  av¬ 
erage  value  0  «  1.8,  as  follows  from  (3-3),  we  get  L_  =  1.8  6,  which 
agrees  with  the  values  given  in  the  majority  of  textbooks  and  handbooks 
on  heat  transfer  [9-11]  for  the  equivalent  radii  for  a  radiating  plane- 
parallel  layer.  In  the  literature  cited  there  are  also  plots  of  the 
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hemispherical  blackness  coefficients  e.  =  e  (p,  T,  L_)  for  water  vapor  . 
and  carbon  dioxide. 

The  result  of  the  calculation  of  by  Formulas  (3.2)  and  (3*3) 
differ  insignificantly  (on  the  high  side)  from  the  results  of  the  cal¬ 
culation  which  takes  into  account  the  configuration  of  the  shock  layer.* 
The  latter  can  be  obtained,  for  example,  by  breaking  up  the  spherical 
segment  Into  n  sectors  with  center  at  the  critical  point,  and  replacing 
in  the  i-th  sector  the  boundary  (the  shock  wave)  by  a  circular  arc 
with  average  radius  r^.  It  is  easy  to  verify  that  In  this  case  the  cal¬ 
culation  formula  assumes  the  form 

n 

7,  =  (3- *0 

Here  and  are  the  angles  between  the  ^  axis  and  the  radius 

vectors  bounding  the  _i-th  sector. 

Formula  (3>^)  enables  us  to  take  Into  account  the  distribution  of 
the  parameters  along  the  length  of  the  radiating  segment.  The  maximum 
deviation  of  the  parameters  in  the  segment  from  their  values  in  the 
vicinity  of  the  retardation  point  will  obviously  be  on  the  line  where 
the  shock  wave  crosses  the  plane  tangent  at  the  critical  point  to  the 
body  around  which  the  stream  flows.  It  is  easy  to  show,  by  using  the 
relations  on  the  discontinuity,  that  these  deviations  are  insignificant 
and,  for  example,  for  a  hemispherical  nose  they  turn  out  to  be  small 
quantities  of  the  order  6/R  *  k,  i.e.,  the  parameters  in  the  radiating 
volume  change  little,  so  that  we  can  use  Formulas  (3*2)  and  (3*3)  to 
calculate  qz>. 

In  the  absence  of  absorbing  gas  in  the  boundary  layer,  the  flux 
cjp  calculated  in  this  manner  is  added  to  the  convective  flux.  On  the 
other  hand,  in  the  case  considered  here,  the  flux  qp,  as  follows  from 
(2.8)  and  (3.1) ,  acts  on  the  convective  flux  and  the  additivity  Is  de- 
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stroyed. 

Under  the  assumption*  of  the  present  section.  Expression  ( 3 . 1) , 
using  the  solution  (2.5)  with  P  =  1  and  the  relations 

Pl—P*  r.  "|  ~~  Cl  («•»  —  "H)  +  *1  ’  c,w  ~  —  ym  +fw’  • 

P«P.[* *  =  i«+  (I. e/-Cfce,  +  cM(< *  —  «.> 

can  be  represented  in  the  following  form 

ww.n-n  i‘.+(i.-«j/'i/r  i.+o.-yr  f  „\ 

’*  *,',(/»'  -2/J«.  \l  «w«»* +•*,<*—«!*>  J  / 

*  =  y(«^V  +  (2-».)*«>  (3-5) 

Here  m-^  and  mg  are  the  molecular  weights  of  the  active  gas  and  of 
the  main  stream,  respectively;  G  is  a  proportionality  coefficient  that 
depends  in  the  general  case  on  the  total  pressure  in  the  mixture  and 

on  the  local  value  of  the  temperature.  Assuming 

ronst,  G(p,T)  =  G{p ) 

we  obtain  from  (2.8),  (2.9),  and  (3*5)  after  simple  transformations 


K'~SW-v5v.’  '  '*> 


(3-6) 


•  • 


Here  is  the  convective  heat  flux  calculated  neglecting  the 

radiation  of  the  gas  released  from  the  surface  of  the  body,  and  Cn  Is 

n 

the  number  of  combinations  of  m  taken  n  at  a  time.  The  Integrals  ^ 
are  universal  functions  that  depend  on  the  intensity  of  the  in-blast 
of  the  active  gas;  let  us  present  the  values  of  Iq,  1^,  ...,  1^  for 
three  values  of  fw. 


'•  1  —  1 

1  «-i  |  «-*  |  »->  | 

—  4  | 

I  — * 

5P 

odd 

I  1  I 

10.975 

3.8588 

2.0310 

1.0048 

0.0185 

0.3096 

0.3011 

0.2333 

0.0966 

0.1394 

0.1194 

0.0146 

0.08255 

0.07392 

0.02526 

0-05578 

0.05137 

0.01638 

tP  - 


4.  Formula  .(3* 6) ,  which  was  obtained  under  the  numerous  simplify¬ 
ing  assumptions  of  the  preceding  section,  illustrates  the  qualitative 
influence  of  the  diffusion  of  the  radiating  impurities  on  convective 
heat  exchange.  It  follows  from  (3*6)  that  the  presence  of  radiating 
and  absorbing  impurities  in  the  boundary  layer  can  lead  to  two  oppos¬ 
ing  effects.  The  second  term  of  (3-6)  determines  the  effects  of  deex¬ 
citation  of  the  boundary  layer,  due  to  the  energy  loss  resulting  from 
the  radiation  of  the  active  gas.  This  effect  leads  in  accordance  with 
(3.6)  to  a  reduction  in  the  convective  heat  exchange.  The  third  term 
is  due  to  absorption  by  the  active  gas  of  energy  radiated  by  the  shock 
layer  and  by  the  surface  of  the  body  in  the  stream.  Unlike  the  deexci¬ 
tation  effect,  the  effect  of  reabsorption  leads,  in  accordance  with 
(3*6),  to  an  increase  in  the  convective  fluxes.  It  is  obvious  that 
this  effect  Is  determined  by  the  magnitude  of  the  radiant  flux  ra¬ 
diated  by  the  shock  layer.  However,  the  influence  of  the  radiation  of 
the  shock  layer  is  not  limited  to  the  foregoing  effect.  In  addition, 
the  flow  in  the  boundary  layer  can  also  be  changed  as  a  result  of  a 
change  in  the  gas  parameters  on  its  outer  boundary  due  to  the  deexclta- 
tlon  of  the  shock  layer,  something  not  taken  into  account  in  the  pre¬ 
ceding  sections.  In  order  to  take  this  effect  into  account,  it  is  suf¬ 
ficient  to  substitute  in  the  formulas  obtained  above  the  value  of  the 
enthalpy  on  the  outer  boundary,  calculated  with  account  of  the  deexci¬ 
tation  of  the  shock  layer.  In  this  case  we  can  neglect  the  deexcita¬ 
tion  of  the  boundary  layer  resulting  from  the  radiation  of  the  par¬ 
ticles  of  the  main  stream  (for  example,  air),  by  virtue  of  the  rela¬ 
tively  low  temperatures  in  the  boundary  layer  as  compared  with  the  tem¬ 
pera  ture  of  the  shock  layer  [(T^/Ts)^  «  1],  i.e.,  the  deexcitation  of 
the  boundary  layer  is  due  only  to  the  radiation  of  the  active  gas. 

The  distribution  of  the  enthalpy  between  the  outer  limit  of  the 
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boundary  layer  and  the  shock  wave,  due  to  the  deexcitation  of  the 
shock  layer,  can  be  readily  determined  by  using  the  method  indicated 
in  [1],  namely:  we  put  in  the  equations  of  the  cited  paper  X  =  |i  =  0. 
Then  the  distribution  of  the  velocities  in  the  shock  layer  is  deter¬ 
mined  by  the  solution  of  Reference  [12],  while  the  temperature  profile 
is  given  by  the  equation 

{TwITtr<  I) 

/  Ij  "  ”  239  ■  “J  V.*.  r,  -  (^«1) 

/  -  -  3*)(-J)*+  p*  -  n 


In  the  first  equation  of  (4.1)  we  can  assume  c^  and  a  constant, 
equal  to  some  values  averaged  over  the  thickness  of  the  shock  layer 
[13).  In  this  case,  as  was  indicated  in  [1]  the  temperature  profile  In 
the  shock  layer  Is  determined  In  final  form  by 

+  „  «*JMV  v_  Jtt—  2V2R6 


For  a  hemispherical  small  nose  we  have  6/R  =  k(l  +  V/5E)”^J 
consequently 


and 


(i-  vT5)/2  Ytm 

By  a  method  that  differs  somewhat  from  that  indicated  in  [1], 
this  solution  was  obtained  in  [14).  It  differs  from  the  solution  (4.2) 
only  in  the  fact  that  the  parameter  contains  a  factor  Eg(0.5T^), 
which  is  close  to  unity  and  takes  approximately  into  account  the  reab¬ 
sorption  in  the  shock  layer*  (t^  is  the  optical  thickness  of  the  shock 
layer).  The  solution  obtained  has  a  singularity  at  y  -»•  0,  but  in  the 
case  tinder  consideration  It  can  be  used  to  determine  the  parameters  on 
the  outer  limit  of  the  boundary  layer,  parameters  contained  in  Formulas 
(2.8)  and  (3*6).  This  singularity  was  considered  previously  in  Refer¬ 
ences  [16,  1,  14 ] .  In  connection  with  the  foregoing  singularity,  and 
with  the  need  for  using  a  finite  thickness  of  boundary  layer  when  tak- 
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ing  into  account  the  deexcitation  of  the  shock  layer,  it  must  be  noted 
that  in  the  vicinity  of  the  critical  point,  by  virtue  of  the  large  vor- 
ticity  of  the  stream  in  the  nonviscous  part  of  the  shock  layer,  the 
solution  of  the  radiating  boundary  layer  should  go  over  into  the  solu¬ 
tion  for  a  radiating  shock  layer,  determined  for  the  same  value  of  the 
stream  line.  At  the  same  time,  as  in  flat  bodies,  the  parameters  on 
the  outer  limit  of  the  boundary  layer  are  equated  to  the  parameters  on 
the  zero  stream  line  of  the  nonviscous  flow,  obtained  without  account 
of  the  boundary  layer  [4],  which  makes  it  possible  in  many  cases  of 
flow  around  flat  bodies  to  carry  out  an  asymptotic  transition  from  the 
solution  for  the  radiating  boundary  layer  into  the  solution  of  the 
problem  for  the  nonviscous  radiating  shock  layer  [18]. 

It  is  obvious  that  in  calculating  the  radiant  flux  qr  from  a  shock 
layer  with  allowance  for  its  deexcitation,  it  becomes  necessary  to  take 
into  account  the  variation  of  the  parameters  over  the  thickness  of  the 
radiating  segment.  In  this  case  qr  can  be  calculated  by  breaking  up 
the  ^-th  sector  in  (3.4)  into  m^  concentric  layers,  summing  then  the 
radiant  flux  from  all  m^  elements  with  allowance  for  their  attenuation 
in  the  preceding  layers.  For  the  case  of  gray  radiation,  the  calcula¬ 
tion  formula  can  be  represented  in  this  case  in  the  form 

*•  mt  i- 1 

9r=  S  n  (*—  eji)(siB*e,—  sin*  «7y  (4.3) 

*.<=i  t— 

Here  =  ee(rij  “  ri  j_i>  pij»  Tij)  is  the  effectlve  blackness 
coefficient  of  a  hemispherical  gas  volume  with  radius  equal  to  the 
thickness  of  the  ij-th  element,  r^j  is  the  radius  vector  of  the  ij-th 
element,  while  and  p^  are  the  temperature  and  pressure  in  the 
ij-th  element. 

Thus,  the  suimnary  effect  ■  connected  with  the  presence  of  the  radia¬ 
tion  field  in  the  gas  stream,  is  brought  about  by  the  deexcitation  of 
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the  radiating  Impurities  of  the  active  gas  In  the  boundary  layer,  lead¬ 
ing  to  a  decrease  in  the  convective  streams,  by  the  absorption  of  ra¬ 
diant  flux  (Jj,  from  the  shock  layer  by  the  impurities,  leading  to  an 
Increase  in  the  convective  streams,  and  by  the  deexcltatlon  of  the 
shock  layer,  which  leads  to  a  reduction  In  ^  and  In  the  temperature 
on  the  outer  limit  of  the  boundary  layer.  In  addition,  the  active  gas 
screens  somewhat  the  surface  of  the  body  against  the  radiant  flux  q^,, 
i.e.,  the  fraction  of  the  energy  absorbed  by  the  gas  Is  not  offset, 
owing  to  the  relatively  low  temperatures  In  the  boundary  layer,  by  its 
radiation  (it  is  obvious  that  Formula  (4.3)  takes  Into  account  this 
screening  effect). 

Consequently,  the  absorption  of  the  active  gas,  which  leads  (by 
deformation  of  the  temperature  profile  in  the  boundary  layer)  to  an 
Increase  In  the  convective  flux,  reduces  somewhat  the  radiant  flux.  On 
the  other  hand  the  deexcitation,  without  causing  an  appreciable  In¬ 
crease  In  the  radiant  flux  to  the  body  qr,  can  noticeably  attenuate 
the  profile  of  the  temperature  In  the  boundary  layer  and  accordingly 
decrease  the  convective  flux. 

Depending  on  which  of  the  foregoing  effects  plays  the  predominant 
role,  the  action  of  the  radiation  field  on  the  temperature  field  in 
the  perturbed  region  can  lead  to  either  an  Increase  or  to  a  decrease 
in  the  summary  flux  to  the  shell  of  the  body  In  the  stream: 

(4.4) 

where  Is  the  absorbing  ability  of  the  surface  of  the  body  (It  Is 
frequently  assumed  that  =  ew,  in  accordance  with  Klrchhoff's  In¬ 
tegral  law),  and  z  is  a  factor  that  takes  Into  account  the  Influence 
of  the  radiation  field  on  the  convective  heat  exchange. 

In  the  case  under  consideration,  the  factor  z  Is  determined  by 
Formulas  (2.8)  or,  with  a  smaller  degree  of  accuracy,  by  (3-6) .  These 
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formulas  obviously  can  be  used  not  only  for  a  qualitative  analysis  of 
the  influence  of  the  diffusion  of  the  radiating  impurities  on  the  heat 
exchange,  but  also  for  numerical  calculations. 

Indeed,  the  representation  of  qw  in  the  form  q^z  Is  quite  con¬ 
venient,  and  if  qWQ  is  determined  with  allowance  for  the  main  proper¬ 
ties  of  the  real  gas,  then  the  integral  quantities  z  can  be  calculated 
with  sufficient  accuracy  under  the  assumptions  made  in  the  second  and 
even  in  the  preceding  sections.  Obviously,  the  variability  of  l ,  the 
deviation  of  the  Lewis  number  from  unity,  or  an  exact  account  of  the 
terms  containing  the  pressure  gradient,  etc.  can  introduce  only  small 
corrections  into  z,  and  cannot  change  appreciably  the  value  of  the 
thermal  flux,  all  the  more  since  all  these  factors  are  taken  into  ac¬ 
count  by  the  first  (main)  factor  of  the  formulas  obtained. 

From  these  formulas  we  can  readily  establish  also  the  influence 

of  the  blackness  coefficient  of  the  surface  e  on  the  increment  of  the 

w 

convective  flux  due  to  the  radiation.  The  decrease  in  e  as  a  function 

w 

of  the  temperature  factor  Tw/Ts  can  lead  both  to  a  decrease  or'  to  an 
increase  in  the  convective  heat  fluxes. 

If  we  neglect  reabsorption  for  gray  radiation  with  VTS  <  5 

the  decrease  in  ew,  as  follows  from  (3-6)  and  (3*3),  leads  to  an  in¬ 
crease  in  the  convective  fluxes,  while  in  the  case  of  T^Tg  >  it 

leads  to  their  decrease.  The  meaning  of  these  inequalities  is  quite 
obvious.  For  example,  when  T «  0  the  first  inequality  is  always  satis- 
fled,  and  leads  to  an  increase  in  the  fluxes  with  decreasing  e  .  In 
this  case  when  =  0  the  increment  of  the  flux  is  approximately 
doubled  as  compared  with  ew  =  1.  The  latter  is  due  to  the  fact  that  in 
the  former  case  the  boundary  layer  is  traversed  twice  by  the  radiant 
flux  ( incident  from  the  shock  layer  and  reflected  from  the  surface  of 
the  body),  whereas  when  ew  =  1  the  radiant  flux  passing  through  the 
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boundary  layer  Is  completely  absorbed  by  the  wall. 

It  Is  obvious  that  with  Increasing  T  and  with  decreasing  e  ,  the 
principal  role  will  be  assumed  by  the  absorption  of  the  thermal  radia- 
tlon  of  the  body  by  the  active  gas,  and  consequently  If  the  second  of 
the  foregoing  Inequalities  Is  satisfied  an  Increase  In  ew  leads  to  an 
Increase  In  the  convective  thermal  fluxes. 
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[Footnotes  ] 


We  note  that  the  use  of  a  more  general  relationship  does  not 
add  additional  difficulties  to  the  solution. 


52  nils  difference  decreases  with  Increasing  optical  thickness 

of  the  shock  layer  and  consequently  with  Increasing  q^. 

55  According  to  [6],  if  reabsorption  Is  taken  Into  account,  the 

right  half  of  (4.1)  assumes  In  the  first  approximation  the 

2i  h 

form  —pact T  and  only  when  «  0  does  it  go  over  Into  — 2aaT  , 

since  0(0)  =  2,  where  0  is  determined  from  (3*2).  It  follows 
therefore  that  for  an  approximate  account  of  reabsorption  it 
Is  sufficient  to  multiply  the  parameter  in  the  solution 

(4.2)  by  p/2.  It  Is  obvious  that  this  factor  has  the  same 
meaning  as  £2(0.5  t^),  and  when  «  0  It  also  turns  Into 

unity.  We  note  that  both  methods  of  taking  reabsorption  Into 
account  can  be  used  only  for  relatively  small  t^.  Large  op¬ 
tical  thicknesses  are  considered  In  (15-17). 
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EFFECT  OF  RATE  OF  HEAT  RELEASE  IN  ELECTRIC  EXPLOSIONS  IN  WATER 
ON  THE  DISTRIBUTION  OF  THE  EXPLOSION  ENEROT 


M.I.  Vorotnikova 
(Novosibirsk) 


In  the  case  of  explosions  In  water  it  can  be  assumed  that  the  en¬ 
ergy  Is  essentially  divided  between  the  energy  carried  away  by  the 
shock  wave  and  the  energy  of  radial  motion  of  the  layers  of  the  medium 
adjacent  to  the  gas  bubble. 

We  Investigate  In  this  work  the  energy  distribution  of  an  under¬ 
water  electric  explosion  among  the  above-mentioned  types  of  motion,  as 
a  function  of  the  Initial  velocity  of  energy  release.  The  initial  rate 
of  energy  release  was  varied  by  connecting  different  Inductances  In 
the  discharge  circuit. 


Fig.  1 


Fig.  2 


An  attempt  to  determine  the  ratio  of  the  energy  of  the  positive 
phase  of  the  acoustic  wave  in  water  in  the  electric  explosion  to  the 
total  energy  stored  In  the  capacitor,  for  different  parameters  of  the 
discharge  circuits,  was  made  in  [1].  Hie  dependence  of  the  radiated 
sound  energy  on  the  total  discharge  energy,  the  Inductance,  and  the 
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resistance  is  described  in  the  article 

[21. 

Description  of  the  setup.  The  dia¬ 
gram  of  the  setup  is  shown  in  Pig.  1. 

f  s  n  is  A  capacitor  bank  1  with  total  capacitance 

A 

Pig.  3.  A)  psec.  1*35  \l£,  charged  with  rectifier  2  to  a 

voltage  of  20  kv,  was  discharged  through  a  spark  gap  3  20  mm  in  length 

in  a  bath  4  at  a  depth  of  140  mm.  The  bath  measured  500  x  500  x  1000 
■a 

mnr  and  was  filled  with  tap  water.  The  inductance  of  the  circuit  was 
varied  between  2  and  32  tih.  The  time  dependence  of  the  current  and 
voltage  on  the  spark  gap  was  recorded  with  an  0K-25  oscillograph.  Ihe 
current  transducer  was  a  noninductive  shunt  5  made  of  nichrome  ribbon 
with  divider;  the  voltage  transducer  was  divider  6,  made  up  of  type  VS 
resistors  without  threading. 

The  pressure  wave  from  the  spark  was  registered  with  an  E0-58M 
oscillograph.  The  pressure  transducer  7  was  a  cylinder  7  ram  in  diameter 
and  10  mm  high,  made  of  polarized  barium  tltanate  ceramic,  soldered 
with  Wood's  alloy  to  a  long  zinc  rod  of  the  same  diameter.  This  system 
was  placed  in  an  aluminum  tube  of  inside  diameter  10  mm  and  the  gap 
filled  with  wax.  The  pressure  transducer  was  calibrated  in  a  shock 
tube.  The  time  resolution  was  on  the  order  of  5  |xsec. 

Simultaneously  with  registering  the  current,  the  voltage,  and  the 
pressure,  we  photographed  the  first  pulsation  of  the  gas  bubble  on  the 
rotating  disk  of  the  photoregister.  The  illumination  was  by  means  of 
an  XPP-1500  gas-discharge  flash  lamp  through  a  window  in  the  bath. 

Dependence  of  the  energy  of  the  positive  phase  of  the  pressure 
wave,  on  the  rate  of  energy  release.  Figure  2  shows  a  time  sweep  of  the 
current  and  the  voltage  of  the  spark  gap.  The  oscillograph  was  trig¬ 
gered  at  the  instant  when  the  air  gap  8  was  closed  (Fig.  1),  when  the 
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voltage  applied. to  the  spark  gap  In  the  water  Is  20  kv.  However,  the 
breakdown  of  the  spark  gap  In  the  water  is  not  Instantaneous,  but  af¬ 
ter  a  certain  delay,  during  which  the  spark  channel  Is  formed.  The  de¬ 
lay  time  varies  and  reaches  10-15  io.se  c.  Hie  discharge  has  the  charac¬ 
ter  of  damped  harmonic  oscillations;  the  period  and  the  damping  decre¬ 
ment  depend  on  the  parameters  of  the  discharge  circuit. 


3/gt  A  /o  The  current  and  voltage  oscillograms 


wat/c 

* 

y 

7 1 

make  It  possible  to  calculate  the  follow- 

r 

. 

ing: 

✓1 

,  |  1.  The  total  energy  released  In  the 

tin"  B 

4* 

■  \ _ i  \  .. 

a  JW  71F  spark  gap,  given  by 


e 

Fig.  4.  A)  W,  Joule/sec; 
B)  1/L,  IT1. 


Experiment  has  shown  that  the  value 
of  the  integral  varies  little  with  the  In¬ 
ductance  In  the  circuit  and  Is  approxl- 
Flg.  5  mately  equal  to  CU^/2,  i.e.,  all  the  en¬ 

ergy  stored  Initially  In  the  capacitor.  In 
the  experiments  W q  =  270  Joules. 

2.  Using  the  oscillograms,  we  can 
plot  the  dependence  of  W(t)  In  Joules  of 
released  energy  on  the  time,  as  shown  In 
Fig.  3,  where  different  values  of  the  In- 

'  * 

Fig.  6.  A)  W,  Joule/sec.  ductance  L  of  the  discharge  circuit  are 
Indicated  on  the  curves.  From  the  Initial  slope  we  determine  the  Ini¬ 
tial  rate  of  energy  release  W. 

For  each  of  the  seven  values  of  the  Inductance  we  made  20  experi¬ 
ments,  two  with  average  delay  time  were  In  error. 

Experiment  shows  that  the  Initial  rate  of  energy  release  was  In- 


versely  proportional  (Fig.  4)  to  the  inductance  of  the  discharge  cir¬ 
cuit,  L. 

The  energy  Wp  of  the  positive  phase  of  the  pressure  wave  is  an  im¬ 
portant  characteristic  of  the  underwater  explosion  and  serves  as  a 
measure  of  the  useful  work  which  can  be  performed  by  the  shock  wave. 
Assuming  spherical  symmetry,  the  work  W(R)  on  a  surface  of  radius  R, 
performed  after  the  arrival  at  the  wave  front  at  the  instant  t  =  0  is 
determined  by  the  expression 

IF  («)  =  4.1  \  rPVIt 

v 

0 

where  U  is  the  stream  velocity  on  a  surface  of  radius  R  and  TJdt  Is  the 
displacement  of  the  surface  during  the  time  dt. 

If  it  is  required  to  calculate  on  the  basis  of  the  experimental 
pressure-time  curves  (the  oscillogram  on  Fig.  5)  the  energy  radiated 
in  a  pressure  wave  of  finite  amplitude,  then  the  stream  velocity  can 
be  eliminated  by  means  of  the  expression 

U  =  (p  -  Pq)/Pqco 

As  a  result,  neglecting  the  hydrostatic  pressure  p^,  we  obtain 

t 

,r<fl)  =  5 

Here  S  is  the  area  of  the  sphere  whose  radius  is  equal  to  the  dis¬ 
tance  from  the  spark  gap  to  the  pressure  transducer;  Pqcq  is  the  wave 
resistance  of  the  water;  t  is  the  duration  of  the  positive  phase  of 
the  pressure  wave,  and  £  is  the  peak  pressure  [3]. 

Experiment  has  shown  that  the  duration  of  the  positive  phase  of 
the  pressure  wave,  At+,  does  not  depend  on  the  rate  of  energy  release 
in  the  spark  gap.  The  quantity  At+  is  a  function  of  the  initial  radius 
of  the  gas  bubble,  whose  dimensions  are  determined  by  the  distance  be¬ 
tween  the  electrodes  of  the  spark  gap  [4]. 
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At  given  parameters  of  the  discharge  circuit,  the  energy  of  the 
pressure  wave  at  a  distance  of  20  cm  from  the  source  amounts  to  only 
a  few  per  cent  of  the  energy  initially  stored  in  the  capacitor. 

The  dependence  of  rj  (the  ratio  of  the  energy  of  the  positive  phase 
of  the  pressure  wave  to  the  energy  initially  stored  in  the  capacitor) 
on  the  initial  rate  of  energy  release  in  the  spark  gap  W  is  shown  in 
Fig.  6. 

The  energy  of  the  radial  motion  of  the  liquid  in  the  gas  bubble 
can  be  expressed  as  the  potential  energy  of  the  expanding  bubble  in 
terms  of  its  maximum  radius  a  and  the  hydrostatic  pressure  pQ  at  the 
explosion  depth  using  the  formula 

V  ■=  1 


The  maximum  radius  of  the  bubble  is  determined  from  the  photo¬ 
graphs,  one  of  which  is  shown  in  Fig.  7  (time  scan  of  the  expansion 
and  compression  of  the  gas  bubble,  arising  during  the  electric  dis¬ 
charge  in  the  water).  The  photograph  was  taken  through  a  slit  by  scan- 


w 


Sfc 


Fig.  7»  mceK  =  msec. 


ning  on  the  rotating  disk  of  the  photoregister.  The  maximum  radius  of 
the  bubble  is  independent,  within  the  limits  of  experimental  error,  of 
the  rate  of  energy  release.  The  result  is  understandable,  since  the 
total  amount  of  energy  does  not  change,  and  the  energy-release  process 
always  remains  sufficiently  rapid  compared  with  the  expansion  of  the 
gas  cavity.  The  energy  of  radial  motion  of  the  liquid  layers  adjacent 
to  the  gas  bubble  amount  to  about  22$  of  the  energy  initially  stored 
in  the  capacitor,  and  is  independent  of  the  initial  rate  of  energy  re- 
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lease . 


In  conclusion,  the  author  is  grateful  to  R.I.  Soloukhin  for  in¬ 
terest  in  the  work. 
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